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Abstract. In this paper, we study the problem of the nonlinear interaction of impulsive 
gravitational waves for the Einstein vacuum equations. The problem is studied in the con- 
text of a characteristic initial value problem with data given on two null hypersurfaces and 
containing curvature delta singularities. We establish an existence and uniqueness result for 
the spacetime arising from such data and show that the resulting spacetime represents the 
interaction of two impulsive gravitational waves germinating from the initial singularities. 
In the spacetime, the curvature delta singularities propagate along 3-dimensional null hyper- 
surfaces intersecting to the future of the data. To the past of the intersection, the spacetime 
can be thought of as containing two independent, non-interacting impulsive gravitational 
waves and the intersection represents the first instance of their nonlinear interaction. Our 
analysis extends to the region past their first interaction and shows that the spacetime still 
remains smooth away from the continuing propagating individual waves. The construction 
of these spacetimes are motivated in part by the celebrated explicit solutions of Khan- 
O Penrose and Szekeres. The approach of this paper can be applied to an even larger class 

of characteristic data and in particular implies an extension of the theorem on formation of 
Ch trapped surfaces by Christodoulou and Klainerman-Rodnianski, allowing non-trivial data 

on the initial incoming hypersurface. 
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1. Introduction 

1.1. Impulsive Gravitational Waves. In this paper, we study spacetime solutions (A4,g) 
to the vacuum Einstein equations 

R pv = o (l) 

representing a nonlinear interaction of two impulsive gravitational waves. Informally, an 
impulsive gravitational spacetime is a vacuum spacetime which contains a null hypersurface 
supporting a curvature delta singularity. Explicit solutions with such properties have been 
constructed by Penrose [31], and its origin can be traced back to the cylindrical waves of 
Einstein- Rosen [10] and the plane waves of Brinkmann [7J . 

Impulsive gravitational waves have been first studied within the class of pp-waves that 
was discovered by Brinkmann [7J, for which the metric takes the form 

g = 2dudr + H(u, X, Y)du 2 + dX 2 + dY 2 , 

and ([I]) implies that 

— — -0 

dX 2 + W 2 ' ^ 
These include the special case of sandwich waves, where H is compactly supported in u. 
Originally, impulsive gravitational waves have been thought of as a limiting case of the pp- 
wave with the function H admitting a delta singularity in the variable u. Precisely, explicit 
impulsive gravitational spacetimes were discovered and studied by Penrose [21] who gave the 
metric in the following double null coordinate form: 

g = -Idudu + (1 - uQ(u))dx 2 + (1 + uQ(u))dy 2 , (3) 

where O is the Heaviside step function. In the Brinkmann coordinate system, the metric has 
the pp-wave form and an obvious delta singularity: 

g = -2dudr - 5(u)(X 2 - Y 2 )du 2 + dX 2 + dY 2 , (4) 
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where S(u) is the Dirac delta. Despite the presence of the delta singularity for the metric in 
the Brinkmann coordinate system, the corresponding spacetime is Lipschitz and it is only 
the Riemann curvature tensor (specifically, the only non-trivial a component of it) that has 
a delta function supported on the plane null hypersurface {u = 0}. This spacetime turns 
out to possess remarkable global geometric properties [30] • In particular, it exhibits strong 
focusing properties and is an example of a non-globally hyperbolic spacetime. 

In a previous paper, we initiated a comprehensive study of impulsive gravitational space- 
times in the context of the characteristic initial value problem. We were able to construct 
a large class of spacetimes which can be thought of as representing impulsive gravitational 
waves parametrized by the data given on an outgoing and an incoming hypersurface such 
that the curvature on the outgoing hypersurface has a delta singularity supported on a 2- 
dimensional slice. Our construction in particular provides the first instance of an impulsive 
gravitational wave of compact extent and does not require any symmetry assumptions. 

1.2. Collision of Impulsive Gravitational Waves. Returning to the explicit examples, 
one of the interesting features of plane gravitational waves is that they enjoy a principle 
of linear superposition provided that the direction and polarization of the waves are fixed. 
This is not the case when one tries to combine two plane gravitational waves propagating 
in different directions. Nonetheless, explicit solutions to the vacuum Einstein equations 
modelling the interaction of two plane sandwich waves have been constructed by Szekeres 
[55] . Khan- Penrose [T7] later discovered an explicit solution representing the collision of two 
plane impulsive gravitational waves. Further analysis of Khan-Penrose solution was carried 
out by Szekeres [%D] . 




Figure 1. The Khan-Penrose Solution 

The Khan-Penrose solution can be represented by Figure 1. The null hypersurfaces {u = 
0} and {u = 0} have delta singularities in the Riemann curvature tensor. In region I, where 
u < and u < 0, the metric is flat and takes the form 

g = —Idudu + dx 2 + dy 2 . 

In region II, where u < and u > 0, the metric is also flat and takes the form 

g = —Idudu + (1 — u)dx 2 + (1 + u)dy 2 . 

Across the null hypersurface {u = 0} between regions I and II, the curvature has a delta 
singularity. In fact, when u < 0, the Khan-Penrose solution coincides with the Penrose 
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solution ([3]) of one impulsive gravitational wave. The region III, where u > and u < 0, is 
symmetric to region II, and the metric takes the form 

g = —Idudu + (1 — u)dx 2 + (1 + u)dy 2 . 

The intersection of the null hypersurfaces u = and u = represents the interaction of the 
impulsive gravitational waves. Thus region IV, where u > and u > 0, is interpreted as the 
region after the interaction. Here, the metric takes the form 

2(1 -u 2 -u 2 yi , , 

q = . -dudu 

^(l - u 2 )(l - u 2 )(uu + - u 2 )(l - u 2 )) 2 ~ 



9 9 . / 1 — u\/l — u 2 — u\/\ — u 2 , 9 1 + u\/l — u 2 + nV 1 — u 2 , 9 

+ (1 - u 2 - u j v — dx 2 H v — dy 2 

y 1 + m y 1 — u 2 + u v 1 — u 2 1 — uy'l — u 2 — uy 1 — u 2 

Even the spacetime is flat and plane symmetric in regions I, II and III, the curvature is 
nonzero and the plane symmetry is destroyed in region IV, signaling that the two plane 
impulsive gravitational waves have undergone a nonlinear interaction. Nevertheless, the 
metric is smooth when u > 0, u > and u 2 + v 2 < 1. Towards u 2 + u 2 = 1, curvature 
invariants blow up and moreover the metric cannot be extended continuously beyond. 

As seen from Q, the Penrose solution of one impulsive gravitational wave in particular 
belongs to the class of linearly polarized pp-w&ves, which takes the general form 

g = -Idudr - H(u)(cosa(X 2 - Y 2 ) + 2 sin aXY)du 2 + dX 2 + dY 2 . 

The constant a is defined to be the polarization of the wave. Thus the Khan-Penrose solution 
represents the interaction of two linearly polarized impulsive gravitational waves with aligned 
polarization. The Khan-Penrose construction was later generalized by Nutku-Halil |29J who 
wrote down explicit solutions modelling the interaction of two plane impulsive gravitational 
waves with non-aligned polarization. These spacetimes have the same singularity structure 
as that of Khan-Penrose. 

Further examples of interacting plane impulsive gravitational waves were constructed via 
solving the characteristic initial value problem with data prescribed on the boundary of 
region IV. This was undertaken by Szekeres [3D] and Yurtsever [43] for the case of aligned 
polarization via the Riemann method. The general case of non-aligned polarization has 
been studied in a series of papers of Hauser- Ernst [3] , [T5] , [TB] by reducing it to the matrix 
homogeneous Hilbert problem. The construction of even more general plane distributional 
solutions for the vacuum Einstein equations that include colliding impulsive gravitational 
waves was carried out in [22], [23] . 

We refer the readers to [12] , [13] , [1] , [6] and the references therein for further description 
and more examples of spacetimes with colliding impulsive gravitational waves. 

The solutions of Khan-Penrose, Szekeres and Nutku-Halil as well as the Hauser-Ernst so- 
lutions are all constructed within the class of plane symmetry. This imposes the assumptions 
that the wavefronts are flat and that the waves are of infinite extent. It has been speculated 
that the singular structure of the Khan-Penrose solution is an artifact of plane symmetry 
[4TJ . Concerning the assumption of plane wavefronts, Szekeres [10] wrote 

The eventual singular behavior is just another aspect of Penrose's result that 
plane gravitational waves act as a perfect astigmatic lens. It is certainly 
false for waves with curved fronts, but such waves may still act as imperfect 
lenses providing a certain degree of focusing and amplification for each other... 
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Clearly a better understanding of the interaction of gravitational waves with 
more realistic wavefronts is a problem of considerable importance. 

A partial remedy has been suggested by Yurtsever jl2] , who did a heuristic study of "almost 
plane waves" and their interactions, allowing waves of large but finite extent. Our present 
paper considers the interaction of impulsive gravitational waves with finite extent and with 
wavefronts having arbitrary curvature. Locally, this in particular includes the case that the 
wavefronts are flat. Nevertheless, even in this special case, we do not require either of the 
waves to be linearly polarized. 

1.3. Interaction of Coherent Structures. The nonlinear interaction of gravitational 
waves in general relativity can be viewed in the wider context of nonlinear interaction of 
coherent structures such as solitons, vortices, etc. in evolutionary gauge theories, nonlinear 
wave and dispersive equations. The completely integrable models KdV [llj, 1-dimensional 
cubic Schrodinger equation jH] and Sine-Gordon equation [1] not only admit individual 
solitary waves, but also exact solutions representing their superposition. In the past, these 
solutions have an asymptotic form of individual propagating solitary waves. For the period 
after nonlinear interaction, which can be described explicitly and typically results in a phase 
shift, a new superposition of new individual propagating solitary waves emerges in the dis- 
tant future. These solutions are analyzed by means of the inverse scattering method. For 
the non-integrable models, our knowledge is much more limited and only partial results are 
available. In those cases, most of the results concerned perturbative interaction of coherent 
structures in the regimes which are either close to integrable or corresponding to interactions 
with high relative velocity or in which one of the objects is significantly larger than the other 
one. In this context, we should mention the work of Stuart on the dynamics of abelian Higgs 
vortices [33j and the Yang-Mills-Higgs equation |34J and the recent breakthrough work of 
Martel-Merle on the nonlinear solitary interaction for the generalized KdV equation [27] , 
[26]. 

Returning to the present work, one of the main challenges in treating the interaction of 
impulsive gravitational waves is their singular nature, i.e., not only do we want to describe 
precisely how gravitational waves affect each other during the interaction, but we also need 
to content with the fact that each impulsive gravitational wave separately is a singular 
object. We should note that partially because of this challenge, no results of this kind are 
available even for semilinear, let along quasilinear, model problems. On the other hand, 
model problems may not be even suitable for studying the phenomena discovered in this 
work since it is precisely the special structure of the Einstein equations that plays a crucial 
role in our analysis and its conclusions. 

1.4. Previous Work on Impulsive Gravitational Spacetimes. In a previous paper 
[23], we studied the (characteristic) initial value problem for spacetimes representing a single 
propagating impulsive gravitational wave. Corresponding to such spacetimes, we considered 
data that have a curvature delta singularity supported on an embedded 2-sphere So iMs on 
an outgoing null hypersurface, and is smooth on an incoming null hypersurface. We showed 
that such data give rise to a unique impulsive gravitational spacetime satisfying the vacuum 
Einstein equations. Moreover, the curvature has a delta singularity supported on a null 
hypersurface emanating from the initial singularity on Sq^ and the spacetime metric remains 
smooth away from this null hypersurface (see Figure 2). 
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Figure 2. Propagation of One Impulsive Gravitational Wave 

1.5. Description of Results in this Paper. In this paper, we begin the study of the 
(characteristic) initial value problem for spacetimes which represent the nonlinear interaction 
of two impulsive gravitational waves. For such a problem, the initial data have delta function 
singularities supported on embedded 2-spheres So,u and S Uaj o on the initial null hypersurfaces 
Ho and H_ respectively see Figure 3). According to the results that were obtained in [25J, 
before the interaction of the two impulsive gravitational waves, i.e., for u < u s or u < u s , 
a unique solution to the vacuum Einstein equations exists, and the singularity is supported 
on the null hypersurfaces emanating from the initial singularities. 




Figure 3. Nonlinear Interaction of Impulsive Gravitational Waves 

Our focus here will be to understand the spacetime "beyond" the first interaction (region 
IV in Figure 3). We will show that the resulting spacetime will be a solution to the vacuum 
Einstein equations with delta function singularities in the curvature on the corresponding null 
hypersurfaces germinating from the initial singularities. Surprisingly, the spacetime remains 
smooth locally in region IV after the interaction of the impulsive gravitational waves. Our 
main result for the collision of impulsive gravitational waves is described by the following 
theorem: 

Theorem 1. Suppose the following hold for the initial data set: 

• The data on H are smooth except across a two sphere So )Ms , where the traceless part 
of the second fundamental form of H has a jump discontinuity. 

• The data on H_ are smooth except across a two sphere S Us fi, where the traceless part 
of the second fundamental form of H_ Q has a jump discontinuity. 

Then 

(a) For such initial data and e sufficiently small, there exists a unique spacetime (Ai,g) en- 
dowed with a double null foliation u, u that solves the characteristic initial value problem 
for the vacuum Einstein equations in the region < u < u*, < u < u^, whenever 
< e or < e. 
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(b) Let I£ u (resp. H Us ) be the incoming (resp. outgoing) null hypersurface emanating from 
So,u s (resp. S U3i o). Then the curvature components cvab = R(za-> e 4> e s> 64) and a AB = 
R(ba, e3, e_B, 63) are measures with singular atoms supported on H_ u and H Us respectively. 

(c) All other components of the curvature tensor can be defined in L 2 . Moreover, the solution 
is smooth away from H_ u U H Ug . 



Our approach relies on an extension of the renormalized energy estimates introduced in 
[25] . As in [25], our concern is not just the existence of weak solutions admitting two colliding 
impulsive gravitational waves, but also their uniqueness. The uniqueness property follows 
from the a priori estimates developed in this paper and leads to strong solutions of the 
vacuum Einstein equations. 

Parts (b) and (c) of Theorem [I] can be interpreted as results on the propagation of sin- 
gularity for a solution that is conormal with respect to a pair of transversally intersecting 
characteristic hypersurfaces. Similar problems have been studied for general hyperbolic 
equations with a much weaker singularity such that classical well-posedness theorems can be 
applied [5], [2]. In the case of second order equations, it is known that no new singularities 
appear after the interaction of the weak conormal singularities. In general, however, a third 
order semilinear hyperbolic equation can be constructed so that new singularities form after 
the interaction of two weak conormal singularities [52] ■ In this paper, we address stronger 
conormal singularities such that in general, even for semilinear hyperbolic systems, only the 
local propagation of one conormal singularity has been proved [28]. For conormal singulari- 
ties of this strength, no general theorem is known to address the interaction of propagating 
singularities even for semilinear, let alone quasilinear, equations. By contrast, in this work, 
the special structure of the Einstein equations in the double null foliation gauge has been 
heavily exploited to show that even for the stronger conormal singularities that we consider, 
the spacetime remains smooth after their interaction. 

In this paper, as in [2H], we prove a more general theorem on the existence and uniqueness 
of solutions to the vacuum Einstein equations that in particular implies Theorem pTa). In 
addition to allowing non-regular characteristic initial data on both Hq and Hq, our main 
existence theorem extends the results in [25] in two other ways. First, we consider the 
characteristic initial value problem with initial data such that the traceless parts of the null 
second fundamental forms and their angular derivatives are only in L 2 in the null directions 
as opposed to being in L°° in the previous work. Second, in [25J, the constructed spacetime 
lies in the range of the double null coordinates corresponding to {0 < u < e} fl {0 < u < e}. 
In this paper, using some ideas in [21], we extend the domain of existence and uniqueness 
to a region that is not symmetric in u and u, i.e., in ({0 <u<e}n{0<u< Ii}) U ({0 < 
u < e} n {0 < u < I2}), where I\ and I2 are finite but otherwise arbitrarily large (see Figure 
4). We refer the readers to Sections 1.7 and [3] for precise formulations of the existence and 
uniqueness theorem. 

One of the unexpected consequences of our approach in this paper is that we can also apply 
it to the problem on the formation of trapped surfaces. The work of Christodoulou [8] was a 
major breakthrough in solving the problem of the evolutionary formation of a trapped surface 
and this was later extended and simplified in [19], [20]. In all of those works, characteristic 
initial data were prescribed on Hq fl {0 < u < e} and Hq with sufficient conditions for data 
on H n {0 < u < e} formulated in such a way as to guarantee the appearance of a trapped 
surface in the causal future of Hq fl {0 < u < e} and Hq (see Figure 5). 
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Figure 4. Region of Existence 




Figure 5. Formation of Trapped Surface 



The sufficient condition on H fl {0 < u < e} required that certain geometric quantities 
are large with respect to e and thus lead to the problem of constructing a semi-global large 
data solution to the Einstein equations. In all those works, to control the dynamics of the 
Einstein equations, the largeness of geometric quantities associated to H D {0 < u < e} was 
offset by requiring the data on H_ Q to be the trivial Minkowski data. 

Our new approach allows us to eliminate the requirement that the data on H_q have to be 
trivial. It can be replaced by a condition that the data on H_ are merely "not too large" and 
still guarantee the formation of a trapped surface in the causal future of Hq D {0 < u < e} 
and Hjq. We refer the readers to Section [8] for a more precise formulation of the theorem on 
the formation of trapped surfaces. 

1.6. A Toy Model. One of the most challenging aspects of the vacuum Einstein equations 
is its quasilinear and tensorial nature. Nonetheless, it may be instructive to examine a 
related phenomenon in a toy model of a scalar semilinear wave equation satisfying the null 
condition in R 3+1 



(or more generally a system = Q($, $), where $ : IR 3+1 — > W 1 and Q(<3>, $) is a null 
form) with the characteristic initial data 



prescribed on the light cones H = {u := t + r = 0} and H = {u := t — r + 2 = 0} 
respectively and 

h{6) = <f)(u = 0,u = 0,9) 
prescribed on the initial 2-sphere defined by {u = 0, u = 0}. 




(5) 



i<3 



/M) 

g(u,e) 



du4>(u,u = 0,9), 
d u (j)(u = 0,u,9), 
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For this toy model, the analogue of the problem addressed in Theorem [T] is the local 
existence and uniqueness result for ^ in the region {0 < u < e} U {0 < u < 1} for the data 

f = fl + l{«-f>0}/2 

and 

9 = 9i + l{ M -i>o}#2, 

where fi, f 2 , gi, g 2 , h are smooth functions and 1 is the indicator function. For these data, 
duf and d u g have delta singularities supported on the 2-spheres {u = 0} fl {u = |} and 
{ik — 0} PI {u = |} respectively. It turns out that the corresponding solution is smooth 
away from the set {u = J} U {u = ^}, but yet dy4> (resp. d u (()) remains discontinuous across 
{u = §} (resp. {u = §}f| 

Theorem [l] is embedded in a more general local existence and uniqueness result (stated 
precisely in Theorem [2] below). Its analogue for the above toy model is the local existence 
for ([5]) with the data /, g and h only satisfying 

I^VIIl 2 (h (o,£)) < c, 



E 
E 

and 



E 



i<A 

where Q G {x\d X2 — x 2 d Xl ,X2d X3 — x 3 d X2 ,x 3 d Xl — Xidx 3 }. The corresponding solution exists 
in the region {0<n<e}U{0<u<l} and obeys the following estimates: 



sup Y^\pd %( j>\y {Hu) <c', 

sup VH^^IU^j^C, 

i<»<e . 
i<3 

sup ^||fi>|| L2(Hti) + sup ^||fi>|| L 2 ( ^j < C7'. 



0<u<l ~~ 0<n<e . 
— — i<A «<4 



Even though this model hardly reflects the difficulties of the nonlinear structure of the 
vacuum Einstein equations, such local existence, uniqueness and propagation of singularity 
results to our knowledge are not known for this type of equations but follow from the methods 
used in this paper. 

1.7. First Version of the Theorem. Our general approach is based on energy estimates 
and transport equations in the double null foliation gauge. This gauge was used in our 
previous work |25j . The general approach in the double null gauge has been carried out in 
HE], [H] and [IS]. 

The spacetime in question will be foliated by families of outgoing and incoming null 
hypersurfaces H u and H_ u respectively. Their intersection is assumed to be a 2-sphere denoted 
by S Ut u. Define a null frame {e±, e^-, e^, e^}, where e 3 and are null, as indicated in Figure 
6, and e±, e 2 are orthonormal vector fields tangent to the two spheres S u ,u. is tangent to 
H u and e 3 is tangent to H^ u . 



Assuming, of course, that the initial data /2 (resp. 32) is non-zero for u = | (resp. u = \). 
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So,o 

Figure 6. The Null Frame 



Decompose the Riemann curvature tensor with respect to the null frame {ei, e 2 , e 3 , e 4 }: 
a AB = R(e A , e 4 , e B , e 4 ), a AB = R(e A} e 3 , e B , e 3 ), 

Pa = ^R(e A , e 4 , e 3> e 4 ), = ^(e^, e 3 , e 3 , e 4 ), 

p = |-R(e 4 , e 3 , e 4 , e 3 ), ^ = 7 *#(e 4 , e 3> e 4 , e 3 ) 

In the context of the interaction of impulsive gravitational waves, the a and a components 
of curvature can only be understood as measures. In the main theorem below, we do not 
require a and a to even be defined. 

Define also the following Ricci coefficients with respect to the null frame: 

Xab = g{D A e A ,e B ), X AB = g(D A e 3 , e B ), 

Va = --g{D 3 e A , e 4 ), r/ A = -^g(D 4 e A , e 3 ) 
w = -^(-D 4 e 3 ,e 4 ), u = --g(D 3 e A ,e 3 ), 
(a = ^g(D A e A , e 3 ) 

Let x (resp. x) De the traceless part of x (resp. x)- F° r the problem of the interaction of 
impulsive gravitational waves, we prescribe initial data on H (resp. H Q ) such that x (resp. 
X) has a jump discontinuity across So )U (resp. S Ust o) but smooth otherwise. 

As mentioned before, we prove a theorem concerning existence and uniqueness of space- 
times for a larger class of initial data than that for the interacting impulsive gravitational 
waves. The following is the main theorem in this paper on existence and uniqueness of 
solutions to the vacuum Einstein equations. 

Theorem 2. Let 9 A be transported coordinates on the two sphere S^jF] and 7 be the spacetime 
metric restricted to S Ut u. Prescribe data such that Q = 1 [^J Suppose, in every coordinate 
patch on H and H Q , 

det 7 > c, 

Ek^)'^i+ Eiffel sc. 

i<4 i<3 



2.2 



2 see definition in Section 

3 for 2fl~ 2 = —g{L\L[), where L' and i' are denned to be null geodesic vector fields (see Section 2.1). 
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On H , 

£ / \(^yxAB\ 2 du+Y,\(§- d ytrx\<c, 

i<3 J i<3 

and on H_ , 

Then for e sufficiently small, there exists a unique spacetime solution (Ai,g) endowed with 
a double null foliation u, u and solves the characteristic initial value problem for the vacuum 
Einstein equations in ({0 < u < e} n {0 < u < h}) U ({0 < u < e} n {0 < u < I 2 }), 
where I\ and I 2 are finite but arbitrarily large. Associated to the spacetime a double null 
coordinate system (u, u, 8 1 , 6 2 ) exists, relative to which the spacetime is in particular Lipschitz 
and retains higher regularity in the angular directions. 

Due to the symmetry in u and u, it suffices to prove the Theorem in < u < I, < u < e. 
In the sequel, we will focus on the proof in this region. The other case can be treated 
similarly. A more precise formulation of the theorem can be found in Section |3j 

In this paper, local existence and uniqueness is proved under the assumption that the 
spacetime is merely W 1,2 . In terms of differentiability, this is even one derivative weaker 
than the recently resolved L 2 curvature conjecture ([21], [35], [36], [37], [38J ) . Of course the 
W 1 ' 2 assumption refers to the worst possible behavior observed in our data and our result 
heavily relies on the structure of the Einstein equations which allows us to efficiently exploit 
the better behavior of the other components. 

Theorem [2] in particular shows the existence and uniqueness of solutions for the initial data 
of nonlinearly interacting impulsive gravitational waves. An additional argument, based on 
the estimates in the proof of Theorem [2j will be carried out to show the regularity of the 
spacetime with colliding impulsive gravitational waves, i.e., parts (b) and (c) in Theorem [T] 
Theorem 2] also forms the basis for the theorem on the formation of trapped surfaces 
In particular, Theorem extends the existence theorem of Christodoulou 



(Theorem |5_ 

[8] to data that is not necessarily small on H_ while allowing the data to be large on H Q . 
Moreover, the estimates obtained in Theorem [2] show that for a large class of data on H_ Q 
that is not necessarily close to Minkowski space, there exists an open set of initial data on 
H such that a trapped surface is formed in evolution. 

1.8. Strategy of the Proof. Without symmetry assumptions, all known proofs of existence 
and uniqueness of spacetimes satisfying the Einstein equations are based on L 2 -type estimates 
for the curvature tensor and its derivatives or the metric components and their derivatives. 
One of our main challenges in [25] and this paper is that for an impulsive gravitational wave 
the curvature tensor can only be defined as a measure and is not in L 2 . 

Let \l/ denote the curvature components and T denote the Ricci coefficients. In [25] where 
we studied the propagation of one impulsive gravitational wave, the curvature component 
a is non-L 2 -integrable. Nevertheless, we showed that the L 2 -type energy estimates for the 
components of the Riemann curvature tensor 

m 2 + I m 2 < / ^ 2 + / m 2 + I j j vmmdu'du. (6) 

h u Jh„ Jh Jh () Jo Jo Js„, „, 



4 In fact, one of the motivations for formulating Theorem [2] for a finite but arbitrarily long u region is for 
proving Theorem [5J 
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coupled together with the null transport equations for the Ricci coefficients 

v 3 r = * + it, v 4 r = * + IT 

can be renormalized and closed avoiding the singular curvature component a. 

In this paper, we consider spacetimes with two interacting impulsive gravitational waves 
and therefore both curvature components a and a are not L 2 -integrable. We thus need to 
extend the renormalization in |25J and to close the energy estimates circumventing both a 
and a. 

In the remainder of this subsection, we will explain the main ideas for proving a priori 
estimates. Note that since we are working at a very low level of regularity, a priori estimates 
alone do not imply the existence and uniqueness of solutions. An additional argument to 
go from a priori estimates to existence and uniqueness was carried out in [25] in which we 
studied the convergence of a sequence of smooth solutions of the vacuum Einstein equations 
to the non-regular solution. A direct but tedious modification of that argument can be 
carried out in the context of this paper, giving the desired existence and uniqueness result. 
We, however, will be content to prove a priori estimates in this paper and refer the readers 
to [2S] for more details. 

After we explain the ideas for proving the a priori estimates, we will then return to sketch 
the ideas in the proofs of the regularity for colliding impulsive gravitational waves (Theorem 
[TJb),(c)) and the formation of trapped surface. 

1.8.1. Renormalized Energy Estimates. In [25], we introduced the renormalized energy es- 
timates for the vacuum Einstein equations. This allowed us to avoid any information of a 
while deriving the a priori estimates. In this paper, since in addition to an incoming im- 
pulsive gravitational waves there is an outgoing impulsive gravitational wave, both a and a 
are non-L 2 -integrable. We thus need to renormalize the curvature components in a way that 
avoids both a and a. 

To this end, we view the vacuum Einstein equations as a coupled system for the Ricci co- 
efficients T and the curvature components which is traditionally treated by a combination 
of estimates for the transport equations for T coupled with the energy estimates for curva- 
ture. The renormalization used in this paper replaces the full set of curvature components 
\l/ with the new quantities 

# = tf + IT for v = p, P ,a,p, 

\& = otherwise. 

We also replace the full set of transport equations for F with a subset which does not 
involve the prohibited curvature components a, a (or rather, involves only the renormalized 
components Similarly, we consider a subset of Bianchi equations. We then show that 
the reduced system can still be closed by a combination of transport-energy type estimates. 

To illustrate the renormalization, we first prove the energy estimates for /3 on H u and for 
(p, a) on H_ u by considering the following set of Bianchi equations: 

V 4 p = div - l -x ■ a + r# , 
V 4 a = -div *P + ■ a + TV, 

V 3 /3 = Vp + vv + r*, 
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where ^ denotes the regular curvature components. However, the curvature component a 
still appears in the nonlinear terms in these equations. In order to deal with this problem, 
we consider the equations for the renormalized curvature components p = p — \\ ' X an d 
a = a + \x A x instead. Using the equation 

Vax = -a + rr, 

we notice that the equations can be rewritten as 

v 4 p = div p + r# + rvr + rrr, 
v 4 <r = -div *p + r\fr + rvr + rrr, 
v 3 /3 = Vp + v*<t + r^ + rvr + rrr. 

We now have a set of renormalized Bianchi equations that does not contain a. Using these 
equations, we derive the renormalized energy estimate 

^2 + / ^2 < / ^2 + / ^2 + / / / r ^ + rvr ^ + rrr^du'dyf, 

h u Jh^ J h Jh Jo Jo Js u , a , 

in which a does not appear in the error term. 

It turns out that the same renormalization p and a that was used to avoid a also can also 
be applied to circumvent a. For example, a enters as source terms in the following Bianchi 
equations, 

V 3 p = -div p - l -x ■ a + rtf, 
V 3( j = -div *p + ^x ■ a + r*. 

Using the equation 

v 3 x = -« + rr, 

we see that a does not appear in the equations for V3/) and V30". 

As a consequence, we obtain a set of L 2 curvature estimates which do not explicitly 
couple to the singular curvature components a and a. We say explicitly that there is still 
a remaining possibility that the Ricci coefficients T appearing in the nonlinear error for the 
energy estimates may depend on a and a 

1.8.2. Mixed Norm Estimates for the Ricci Coefficients. In order to close the estimates, it 
is necessary to obtain control of the Ricci coefficients via the transport equations 

v 3 r = * + rr, v 4 r = * + rr. (7) 

In we showed that T can be estimated in L°° by considering a subset of the transport 
equations that do not involve the singular curvature component a (and involve only the 
renormalized curvature components ^). 

In the setting of this paper, in addition to proving bounds on T without any information 
on both singular curvature components a and a, an extra challenge is that unlike in [25], not 
all Ricci coefficients are bounded in the initial data. In fact, for the class of initial data con- 
sidered in this paper, x (resp. x) is om y assumed to be in L"^H 3 (S) (resp. L^H 3 ^)), where 
H 3 (S) refers to the L 2 norm of the third angular derivatives on the 2-spheres. Therefore, 
(uj) at best implies that x (resp. x) can be estimated in L^L^°L 00 (S') (resp. L 2 L^ 'L 00 '(S)) , 
where the L°° norms on the sphere and along the u (resp. u) direction are taken first, before 
the L 2 norm in u (resp. u) is taken. 
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Because of the weaker assumption on the Ricci coefficients in the initial data, we only 
prove estimates for the Ricci coefficients in mixed norms. In fact, we prove different mixed 
norm bounds for different Ricci coefficients. Using a schematic notation if> G {tr%, trx, rj, 77}, 
V>h e {xM and ipH E {x,u}, we only control if) in L™L™L°°(S), in L^L™L°°{S) and 
V» S in L*2£L~ (5). 

It is a remarkable fact that the Einstein equations possess a structure such that these 
mixed norm bounds are sufficient to close all the estimates for the Ricci coefficients using 
the transport equations, as well as the energy estimates for the curvature components. 

As an example, in order to estimate if) in L£° L°° (S) , we use the transport equation 

V 3 ip = l + p + W> + (V' + V'if)(V>+V>ff)- 

Notice that the term if>H does not appear as the source of this equation. Therefore, with the 
control of the Ricci coefficients in the mixed norms, all terms on the right hand side can be 
bounded after integrating in the e 3 (i.e. u) direction to obtain the desired bound for if). 

On the other hand, the transport equation for ipjj contains both ipjj and ipn in the inho- 
mogeneous term: 

Integrating this equation, we get 

\\iPh\\l%>l°°(S) < Initial Data + \\if> + iPh\\lil°°(s)\\?P + iPh\\l™l°°{s) + ■■■■ (8) 

The initial data term and the factor H^Hl^l^s) are not bounded. Nevertheless, since we 
are only aiming to prove estimates for if>n in L^L^ L°° (S) , we can take the L 2 U norm in rfsj) 
and every term on the right hand side is controlled by the mixed norms. This allows us to 
prove the mixed norm estimates for all the Ricci coefficients. 

Even more remarkable is that the bounds we obtain for the Ricci coefficients in mixed 
norms are also sufficient to close the energy estimates for the renormalized curvature com- 
ponents. Schematically, the renormalized energy estimates read as follows: 

\\(^,p,a)\\ L ^^ L 2 {s) + \ \(p,cr,f3)\\ L ^ L 2 L 2 {s) 

< Initial Data + \\T^/ 2 \\ L i L ^ L i {s) + ||r 5 || L i^ x i (S ) + ... 

The error terms on the right hand side have to be controlled by the L 2 curvature bounds on 
the left hand side together with the estimates for the Ricci coefficients in the mixed norms. 
As an example, an error term ipHPP can be controlled after applying Cauchy-Schwarz as 
follows: 

W^hPp\\lili_lhs) < \\p\\Lgi%i?{S)\Wn\\i%L°°L°°{S)- 

Here, it is important to note that using the mixed norms for if>H, we can estimate in L°° 
first, before taking the L 2 norm. On the other hand, an error term of the type ipaPfl cannot 
be controlled in L^L^L 1 ^) since each of the three factors can only be bounded after taking 
the L 2 U norm. Miraculously, such terms never arise as error terms in the energy estimates! 
A similar structure also arises in the error terms of the form 

||r 5 || L iLlLi(5)- 

For this term, ifijj (or ipn) appears at most twice, allowing us to estimate each of them in 
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In order to close all the estimates, we need to prove mixed norm estimates for higher 
derivatives of the Ricci coefficients and energy estimates for higher derivatives of the cur- 
vature components. This is achieved using only angular covariant derivatives V as commu- 
tators. For such estimates, the singular curvature components a and a never arise in the 
nonlinear error terms. Moreover, there is a structure similar to that described above for the 
higher order estimates that allows us to close merely with the mixed norm bounds. 

1.8.3. Estimates in an Arbitrarily Long u Interval. In our main theorem, we prove existence, 
uniqueness and a priori estimates in a region such that only the u interval is assumed to be 
short, while the u interval can be arbitrarily long (but finite). This poses an extra challenge 
since when we control the nonlinear error terms integrated over the u interval, we do not 
gain a smallness constant. 

This difficulty already arises in the problem of existence in such a region with smooth 
initial data. This was studied in [HQ It was noticed that both in carrying out the Ricci 
coefficient estimates and the energy estimates for the curvature components, the structure 
of the Einstein equations allows us to prove that whenever a smallness constant is absent, 
the estimate is in fact linear. 

To achieve the bounds of the Ricci coefficients, the following structure of the null structure 
equations was used. Let 

Ti G {trx,X,trx,X,^,^}, r 2 = 77, L 3 = u. 
They satisfy the following transport equations: 

v 4 r! =* + (ri + r 2 + r,)^ + r 2 + r 3 ), 

v 3 r 2 =* + (r x + r 2 )r 1; (9) 

v 3 r 3 =^ + (r 1 + r 2 + r 3 )(r 1 + r 2 ). 

We prove the bounds for T\, T 2 , T 3 in the setting of a bootstrap argument in which 
the control for the curvature components \I> is assumed. The estimates for Ti can easily 
be obtained since integrating in the e 4 (i.e., u) direction gives a smallness constant. For 
T 2 , the integration is in the e 3 (i.e., u) direction and does not have a smallness constant. 
Nevertheless, using the bounds for Ti that have already been obtained, the error term is 
linear in T 2 ! This can thus be dealt with using Gronwall's inequality. Finally, the equation 
for T 3 is also linear in T 3 . Therefore, using the the estimates already derived for Ti and T 2 
together with Gronwall's inequality, the equation for T 3 can be applied to get the desired 
control for T 3 . 

In the energy estimates for the curvature components, there is likewise a term without a 
smallness constant. Nevertheless, it was noted in [21] that the only term not accompanied 
by a smallness constant is also linear. Thus, as in the case in controlling the Ricci coefficient, 
the energy estimates can be closed using Gronwall's inequality. 

Returning to the setting of this paper, this challenge of having an arbitrarily long u interval 
is coupled to the difficulty that the curvature components a and a are singular and that the 
Ricci coefficients x, x, <±L can only be estimated in appropriate mixed norms. As a result, 
unlike in [24] , we cannot use the V4 equations for x an d tr% to gain a smallness constant. 
The V4X equation is unavailable because a appears as the source of this equation, and in 

5 In |24j , the a priori estimates were proved in the case where the u interval is assumed to be short and the 
u interval is allowed to be arbitrarily long. We outline the main ideas of [24] assuming instead the setting in 
this paper. 
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this paper, due to the singularity of a, one of our goals is to prove all estimates without 
any information on a. The V 4 tr\; equation, while can be used, has \x\ 2 as a source term. 
Since \ can only be estimated in L\ using the mixed norm bounds, the integration in the u 
direction does not give a smallness constant. 

Nevertheless, a different structure can be exploited to overcome this challenge. We group 
the Ricci coefficients into Ti, T 2 , T 3 and T 4 according to the equations and estimates that 
they satisfy. Let 

Ti e {trx, x, v, v}, r 2 = 77, r 3 e {%, u}, r 4 = tr%. 

They satisfy the following transport equations: 

v 4 i\ = # + (r\ + r 2 + r 3 + r 4 )(r\ + r 2 + r 3 + r 4 ), 
v 3 r 2 = * + (r 1 + r 2 )r 1> 
v 3 r 3 = * + (i\ + r 2 + r 3 + r 4 )(r x + r 2 ), 
v 4 r 4 = (r 3 + r 4 )(r 3 + r 4 ). 

Notice that T 3 corresponds to the Ricci coefficients ipH and can only be estimated in 
LlL^L^(S). ' 

As before, the control of $ is assumed in a bootstrap setting. The equations for Ti and 
T 2 have similar structures as ([9]). Thus, we first estimate Ti, using the smallness constant 
provided by the integration in the e 4 (i.e., u) direction. We then control T 2 noting that 
with the bounds already obtained for r 1; the error term is linear in T 2 . The equation for 
T 3 is similar to (j9|, except for an extra term containing T 4 , which has not been estimated. 
Nevertheless, T 3 are the terms \ an d oj which are only estimated in L^L™ L°° (S) . Thus the 
error term containing T 4 only has to be controlled after taking the L\ norm. This provides an 
extra smallness constant. Finally, while T 4 satisfies an equation in the e 4 (i.e., u) direction, 
T 3 r 3 appears as a source. Recall that since T 3 can only be controlled in L 2 L^ L°° (S) , this 
error term is only bounded in L l u L^ L°° (S) . In other words, integrating this equation does 
not give a smallness constant. Nevertheless, we can use the control for T 3 derived in the 
previous step! Thus we obtain the desired bounds for all the Ricci coefficients. 

In a similar fashion, the energy estimates also have to be carried out in two steps. Recall 
from ^ that in establishing the energy estimates, we need to control the error terms 

lir**ii W i(s), 

where ^ are the renormalized curvature components. The most difficult error terms are 
those containing (3. This is because /3 can only be controlled in L 2 (H). In order to control 
the error terms, the L 2 (H) norm of (3 has to be integrated over the long w-interval and the 
estimates do not have a smallness constant. To deal with this problem, we first control ft in 
L 2 {W) and (p, a) in L 2 (H). While deriving these bounds, all the error terms are accompanied 
by a smallness constant ea. We estimate /3 after we obtain these bounds. The error terms 
that contain /3 are 

\\xm\LlLlLHsf\ 

and 

WxPPWlilil^s)- 

6 To be more precise, the term that actually appears is ||x/3Vx||liz,ili(s)- We note that using elliptic 
estimates, the control for Vx can be retrieved from the bound for j3. We omit the technical details here. 
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Since the f3 has been controlled first, the first error term is sublinear. For the second term, 
it can be shown that the estimates for x is independent of the bounds on the curvature and 
this term is therefore a linear term. It can thus be dealt with using Gronwall's inequality. 

1.8.4. Signature. In the proof of the a priori estimates, the structure of the Einstein equations 
plays a crucial role. It is thus useful to understand the structure of the equations in a 
more systematic fashion. Here, inspired by the work of Klainerman-Rodnianski [19] on the 
formation of trapped surfaces, we introduce a notion of signature that allows us to explain 
and tract that certain undesirable terms do not appear in a particular equation. Such a 
notion of signature is intimately tied to the scaling properties of the Einstein equations. 

1.8.5. Nonlinear Interaction of Impulsive Gravitational Waves. As mentioned above, The- 
orem [2] implies the existence and uniqueness of solutions to the vacuum Einstein equations 
with characteristic initial data as in Theorem [TJ In the setting of the nonlinear interaction of 
impulsive gravitational waves in Theorem [TJ however, the initial data are more regular than 
the general initial data allowed in the assumptions of Theorem [2] In particular, on each of 
the initial null hypersurfaces, the initial data are only singular on an embedded 2-sphere. 
This allows us to prove that the spacetime is smooth away from the null hypersurfaces em- 
anating from the initial singularities. Moreover, a and a can be defined as measures with 
singular atoms supported on these null hypersurfaces. 

We first note that standard local well-posedness theory and the results of [25] imply that 
the spacetime is smooth in {0 < u < u s } U {0 < u < u s }. Thus in order to show that the 
spacetime is smooth away from the null hypersurfaces {u = u s } and {u = u s }, we only need 
to demonstrate the regularity of the spacetime in {u > u s } fl {u > u s }. 

It turns out that using the a priori estimates derived in the proof of Theorem |2j this can 
be shown by directly integrating the null structure equations. For example, while V 4 x has 
a delta singularity across u = u s , we can prove that it is bounded for u>u s . To this end, 
we consider 

Commute the equation with the V4 derivative and substituting appropriate null structure 
equations, we get 

V3V4X + ^trxV 4 X - 2oA7 4 X 

where ... denotes terms that have already been estimated in the proof the Theorem [3j Thus 
by integrating this equation, we conclude that V4X inherits the regularity of the initial data 
and is bounded as long as u 7^ u s . This procedure can be carried out for all higher derivatives 
to show that the spacetime is smooth in the region {u > u s } fl {u > u s }. 

A surprising feature of this proof of smoothness of the resulting spacetime is that it does 
not require a and a to have delta singularities supported on the corresponding 2-spheres. 
In fact, if the initial data satisfy the assumptions of Theorem [2] and are more regular for 
u > u on Hq and u > u on if , then the spacetime can be proved to be more regular in 
{u > u} H {u > u}\ 

Returning to the interacting impulsive gravitational waves, we show that a and a can be 
defined as measures with delta singularities supported on E_ u and H Us respectively. To see 
this, consider the equations 

a = -V 4 % - tixx ~ 2^X, 
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and 

a. = -V 3 X - trxx - 2u_x- 

We can prove that \ (resp. x) is smooth except across u = u s (resp. u = u s ) where it has 
a jump discontinuity. This implies that a and a are well-defined as measures and they have 
delta singularities supported on H_ u and H Ug respectively. 

1.8.6. Formation of Trapped Surface. Using the existence and uniqueness result in Theorem 
[2j we construct a large class of spacetimes such that the initial data do not contain a trapped 
surface, and a trapped surface is formed in evolution. In particular, unlike in [8], [19] and 
[20], our construction does not require the initial data on I£ Q to be close to that of Minkowski 
space. 

The challenge in this problem lies in the fact that in order to have a trapped surface, certain 
geometric quantities are necessarily large. Recall that in the setting of Christodoulou [8] (see 
Figure 7), characteristic initial data were prescribed on if and a short region of Hq, where 
< u < e. 




Figure 7. Formation of Trapped Surface 
In view of the equation 

V 4 trx = -^(trx) 2 - |x| 2 , 
in order that for some u, trx becomes negative after integrating in a u length of e, x has to be 

1 3 

of size ~ e~2 and consequently a has to be of size ~ e~5. In the work of Christodoulou [H], 
and the later extensions of Klainerman-Rodnianski [19], [20], this largeness of the geometric 
quantities is compensated by requiring smallness of initial data on H_ . 

To go beyond the requirement of Minkowski data on H_ , we notice that while the L^L 00 (S) 
norm of x is large in terms of e, its L^L°°(S') is merely of size ~ 1 with respect to e. 
Therefore, Theorem [2] implies the existence and uniqueness of a spacetime solution for this 
type of initial data, even without any smallness assumptions on H_ Q . Note in particular that 
the assumptions of Theorem [2] do not require any control of a for the initial data. It thus 
remains to show that one can find initial data which do not contain a trapped surface and 
such that a trapped surface is formed in evolution. 

With the initial data that he imposed, Christodoulou identified a mechanism for the 
formation of trapped surface [8]. In view of the equation 

V 4 trx = -^(trx) 2 - |x| 2 , 
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for e sufficiently small, a trapped surface forms in evolution if at u = 0, 

trxO = 0,u = 0,1?) > f \x\ 2 (u = 0,$)du, 

Jo 

while at u = it*, 

trx(u = u*,u = 0,i9) < / |x| 2 ( M = $)ob. 
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To achieve (10) and (11), consider the equations 



and 



V 3 X + ^trxx = V®7? + 2ux - ^trxx + V®V, 



Vstrx + -tr^trx = 2wtrx + 2p — x • X + 2div rj + 2|?y| 



Assuming the right hand side of these equations to be error terms, we get 



V 3 |x| 2 + trxlxl 



0. 



and 



which imply 



and 



V 3 trx + -trxtrx ~ 0, 



|x| 2 (w, u, i?) ~ |x| 2 ( M — 0, u, "&) exp(- 



trx(u, u, $)du') 
trx(u', u — 0, $)du) 



trx(u, u = 0, i?) ~ tix(u = 0, u = 0, i?) exp(— ; 

Christodoulou showed that in the setting of [8], 

trx(w, u, #) ~ trx(-u, w = 0, 

which implies that 

\x\ 2 ( u i "Mii ~ lx| 2 ( u = 0) i?) exp(— / trx(u',u = 0, $)du'). 



(10) 



(12) 
(13) 



(14) 



(15) 



(16) 



Comparing (14) and (16), since tr% < 0, |x| 2 has a larger amplification factor than tr%. 
Therefore, there is an open set of initial data such that a trapped surface is formed in 
evolution. 

In our setting where we remove the smallness assumptions on the data on if , the estimates 
derived in Theorem |4| imply that (12) and (15) hold. Nevertheless, the approximation (13) 



is not necessarily valid. Instead, we impose a condition (50) on H_ Q in Theorem |5| in order 
to guarantee that a trapped surface is formed in evolution. This condition guarantees that 
there is a choice of initial data on H such that ( 10 ) and ( 11 ) hold in the resulting spacetime. 
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Figure 8. The Basic Setup 

2. Setting, Equations and Notations 

Our setting is the characteristic initial value problem with data given on the two charac- 
teristic hypersurfaces H and H_ intersecting at the sphere Sq^q. The spacetime will be a 
solution to the Einstein equations constructed in a neighborhood of H and H_ containing 
So,o- 

2.1. Double Null Foliation. For a spacetime in a neighborhood of S ,o, we define a double 
null foliation as follows: Let u and u be solutions to the eikonal equation 

gTd^vduU = 0, gTdpvdvU = 0, 

satisfying the initial conditions u = on H and u = on H_ . Let 

= -2g^d v u, = -Wd v u. 

These are null and geodesic vector fields. Let 

2n- 2 = -,?(L',L'). 

Define 

e3 = QlJ, e4 = QL' 
to be the normalized null pair such that 

#(63,64) = -2 

and 

L = Q 2 IJ, L = Q 2 L' 

to be the so-called equivariant vector fields. 

In the sequel, we will consider spacetime solutions to the vacuum Einstein equations in 
the gauge such that 

Q = 1, on H and H_ . 

We denote the level sets of u as H u and the level sets of u and H_ u . By virtue of the eikonal 
equations, H u and H_ u are null hypersurface. The sets defined by the intersections of the 
hypersurfaces H u and H_ u are topologically 2-spheres, which we denote by S UjU . Notice that 
the integral flows of L and L respect the foliation S Ut u. 
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2.2. The Coordinate System. On a spacetime in a neighborhood of So,o> we define a 
coordinate system (u, u, 9 1 , 2 ) as follows: On the sphere S^o, define a coordinate system 
(#\# 2 ) such that on each coordinate patch the metric 7 is smooth, bounded and positive 
definite. Then we define the coordinates on the initial hypersurfaces by requiring 

d d 
—9 A = on H Q , and — 9 A = on H . 
ou au 

We now define the coordinate system in the spacetime in a neighborhood of by letting 
u and u to be solutions to the eikonal equations: 

g^d^u = 0, g^d^u = 0, 

and define 6 l ,6 2 by 

fi L A = 0, 

where £l denote the restriction of the Lie derivative to TS Ut u (See [SJ, Chapter 1). Relative 
to the coordinate system (w, u, 9 1 , 9 2 ), the null pair e 3 and e 4 can be expressed as 

for some b A such that b A = on I£ , while the metric g takes the form 

g = -2tt 2 (du ®du + du® du) + r ) AB {dB A - b A du) <g> (d0 B - b B du). 

2.3. Equations. We will recast the Einstein equations as a system for Ricci coefficients and 
curvature components associated to a null frame e 3 , e 4 defined above and an orthonormal 
frame ei,e 2 tangent to the 2-spheres S UiM . Using the indices A, B to denote 1,2, we define 
the Ricci coefficients relative to the null fame: 

Xab = g(D A e 4 ,e B ), x AB = g(D A e 3 , e B ), 

Va = ~-g(D 3 e A , e 4 ), i]_ A = -^g(D 4 e A , e 3 ) 

w = -^g(D 4 e 3 ,e 4 ), u = -^g(D 3 e 4 ,e 3 ), ^ 

(a = ^g(D A e 4 , e 3 ) 
where D A = D e(A) . We also introduce the null curvature components, 

a AB = R(e A , e 4 , e B , e 4 ), a AB = R(e A , e 3 , e B , e 3 ), 

1 1 

/3a = ^R{e A , e 4 , e 3 , e 4 ), (3_ A = -R(e A , e 3 , e 3 , e 4 ), ( lg ) 

p = |i2(e 4 , e 3 , e 4 , e 3 ), er = -jj- *i?(e 4 , e 3 , e 4 , e 3 ) 

Here *i? denotes the Hodge dual of R. We denote by V the induced covariant derivative 
operator on S Uy u and by V 3 , V 4 the projections to S Uj u of the covariant derivatives D 3 , D 4 
(see precise definitions in [IS]). 
Observe that, 

w = --V 4 (logfi), u; = -^V 3 (logft), 

^ = CA + V A (logfi), ^ = -Ca + V^(logfi) 
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Let (j)^ ■ (p^ denote an arbitrary contraction of the tensor product of <jP^ and (f)^ with 
respect to the metric 7. We also define 

(0 (1) §0 (2) )ab := tfV* + <t> { B<t>A - <W0 W ■ (2) ) for one forms 0^ , <f>%\ 



A <p (2) ) A B ■= £ AB (7 _1 ) CD Ac0sd for symmetric two tensors 



where ^ is the volume form associated to the metric 7. For totally symmetric tensors, the 
div and curl operators are defined by the formulas 

(div <f>) Al ...Ar ■= V B 0BAi...A r , 
(curl (f>) Al ...Ar := i 30 ^ ' B<l>CA 1 ...A r . 

Define also the trace to be 

{^4>) Al ...A r ^ ■= (7" 1 ) BC 0BCA 1 ...A r _ 1 . 

We separate the trace and traceless part of x an d X- Let x an d X be the traceless parts 
of x and x respectively. Then x an d X satisfy the following null structure equations: 

V 4 tr X + ^(tr X ) 2 = -|x| 2 - 2wtr X 
V 4 X + trxx = -2ux - a 

V 3 trx+ ^(trx) 2 = -2wtrx- |x| 2 

V 3 X + trx X = - 2 ^X - a 
V 4 trx + ^tr^trx = 2wtrx + 2p - x • X + 2div 77 + 2 1 77 1 2 ( 20 ) 

V 4 X + ^trxx = V% + 2wx - ^trxx + 9% 

V 3 trx + ^tr^trx = 2wtrx + 2p - x • X + 2div r] + 2\r]\ 2 

V 3 X + ^trxx = V®7? + 2wx - ^tr^X + V®V 

The other Ricci coefficients satisfy the following null structure equations: 

V 4 t? = -X-(V-V)-I3 

V 3!Z = -x- (v-v) + P 

V4U = 2uu + ^\v ~ v\ 2 ~ k(v - V) ■ (v + V) - + v\ 2 + \p 
4 4 o 2 

3 1 11 

V 3 w = 2ww + -|?7 -rf + -(77 -77) • (77 + 77) - -\r] + rf + -p 
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The Ricci coefficients also satisfy the following constraint equations 
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div x = ^Vtrx - ^(?7 - V) • (x ~ \^x) ~ P, 

div x = Jvtrx + \{r) -rj)-(x~ o tr x) + P 

- 2 2 2 (22) 

curl i] = —curl r] = a + ~x A x 
K = -p + )pc ■ X - ^trxtrx 



with K the Gauss curvature of the spheres S Ujy ,. The null curvature components satisfy the 
following null Bianchi equations: 



V 3 a + ^trxa = V®/3 + Aua - 3(xP +* X?) + (C + *v)®P, 

V 4 /3 + 2trx/3 = div a - 2w/3 + 77a, 

V 3 /3 + trx/3 = Vp + 2w/3 +*Va + 2 X -/3_ + S(vP +* W), 

V 4 a + *-tT X v = -div */3 + l -x * a - C * P - 2R •* (3, 



(23) 



V 3 a + hrxa = -div + l - X ■* a - ( ■* p - 2 V ■* £, 

V 4 p + ^tr X p = div (3 - l -x ■ a + C • P + 2 R ■ (3, 

3 1 
V 3 p + 2 tr XP = ~ div P ~ 2 X ' - + $ ' & ~ 2r l ' 

V 4 /3 + tr X P = -Vp +* Vd + 2u(3_ + 2x ■ P - S(vp -* rp), 
V 3 /3 + 2tix§_ = -div a - 2w/3 + r?_ • a, 

V 4 a + ^trxa = -V®/3 + 4wa - 3(xp -* x°) + (C - 4r?_)®/3 



where * denotes the Hodge dual on S u> u. 

We now define the renormalized curvature components and rewrite the Bianchi equations 
in terms of them. Let 



P = P~ ljX'% ^ = ^+2^ A i- 
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The Bianchi equations expressed in terms of p and a instead of p and a are as follows: 
V 3 /3 + trx/3 =Vp +* W + 2u(3 + 2 X ■ §_ + 3( V p +* r,&) + ^(V(* • x) -* V(* A *)) 

V 4 cr + ^tr X <x = - div */3 - C •* P - 2 R ■* (3 + l - X * (V%) + •* (5%), 

3 1 ^ 1 ^ 1 

V 4 p + -trxp =div /3 + C- /5 + 2?7-/? — -x- V^t? - • (77(8)77) + -trvjx| 2 , 

V 3 a + ^frycr = - div */3 - C •* §_ - 277 •* §_ + •* (V§77) + ^x •* (^), 

3 1 x ^ 1 

V 3 p + -trxp = - div + ( ■ §_ - 2rj ■ §_ - -% • V®7/ - -% • (778)77) + -trx|x| 2 , 

V 4 ^ + trx^ = - Vp +* V<7 + 2u§_ + 2x ■ P - 3(^/5 -* R a) - ^(V(x • x) -* V(x A x)) 
-3(t/x-X+* VX a x), 

(24) 

Notice that we have obtained a system for the renormalized curvature components in which 
the singular curvature components a and a do not appear. 

In the sequel, we will use capital Latin letters A e {1,2} for indices on the spheres S Ut u 
and Greek letters p G {1, 2, 3, 4} for indices in the whole spacetime. 

2.4. Signature. In this subsection, we introduce the concept of signature. This will allow 
us to easily show that some undesirable terms are absent in various equations. 

To every null curvature component a, (3, p, a, (3, a, null Ricci coefficients x, C, Vi Vi w , and 
the metric components 7, Q, we assign a signature according to the following rule: 

sgn(<f>) = 1 • N 4 (<f>) + (-1) • N 3 (<f>) - 1, 

where iV 4 (0), iV 3 (0) denote the number of times e 4 , respectively e 3 , which appears in the 
definition of <fi. Thus, 

sgn(/3) = 1, sgn(p,a) = 0, sgn(/3) = -1. 

Also, 

sgn(x) = sgn(u) = 1, sgn((, rj : 77) = sgnfr, Q) — 0, s^n(x) = s^n(w) = -1. 

We use the notation and to denote the renormalized curvature component and 
Ricci coefficient respectively with signature s. Then all the equations conserve signature in 
the following sense: The null structure equations are all in the form 

v 4 r^ = \i>( s+1 ) + ^ r( si ) • r^ S2 \ 

Sl+S 2 =S+l 

v 3 r (s) = ^(■'-v + ^2 r (si) • r (s2) - 

«1+S2=S— 1 

and the null Bianchi equations are of the form 

v 4 \t> (s) = v^ (s+1) + ^2 (r (si) • ^ (S2) + r (si) • vr (S2) ), 

Sl+S2=S + l 
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v 3 \i/ (s) = v\& (s_1) + (r (si) • \i> (S2) + r (si) ■ vr (S2) ). 

Sl+S 2 =S—l 

2.5. Schematic Notation. We introduce a schematic notation as follow: Let </> denote an 
arbitrary tensorfield. For the Ricci coefficients, we use the notation 

ip E {trx,trx,^,^}, ipH e ipH e (25) 

Notice that ipH has signature 1 and ipH_ has signature —1. Unless otherwise stated, we will 
not use the schematic notation for the renormalized curvature components but will write 
them explicitly. 

We will simply write ipip (or ipipH, ipP, etc.) to denote arbitrary contractions with respect 
to the metric 7. V will be used to denote an arbitrary angular covariant derivative. The use 
of the schematic notation is reserved for the cases when the precise nature of the contraction 
is not important to the argument. Moreover, when using this schematic notation, we will 
neglect all constant factors. 

We will use brackets to denote terms with any one of the components in the brackets. For 
example, ip(p, a) is used to denote either ipp or ipa. 

The expression VV' 7 will be used to denote angular derivatives of products of Ricci coef- 
ficients. More precisely, V*^ denotes the sum of all terms which are products of j factors, 
with each factor being V %k ip and that the sum of all z^'s being i, i.e., 

^ V v ' 

j factors 



l\+%2 + —+%i 



Using these notations, we write all the equations from Section 2.3 in the schematic form. 



The structure of the equations can be read off directly from Section 2.3 On the other 



hand, we notice that the structure for most of the equations also follows from signature 



considerations as indicated in Section 1.8.4 We will later point out places where we need to 



use an additional structure of the equations that goes beyond signature considerations. 



We first write down the null structure equations (20) and (21) in schematic form. Here, 
we do not write down the two equations that involve the singular curvature components a 
or a. 



V 4 trx = 


Xx + i[)(i/; + il) H ) 




V 3 trx = 


xx + ip(i> + ipH) 




V 4 trx = 


p + Vr] + ip(ip + 4>h) 




V 3 trx = 


p + Vr] + ip(ip + 4>h) 




V 4 r/ = 


P + i/)(i/> + i/>h) 




v 3 r? = 


§_ + (v + v)(tix + 4>H) 


v 4 x = 


p + v?7 + i>{i> + i>u) 


+ i)H_(t?X + 


V 3 x = 


f> + Vp + 1p(lp + 1pn) 


+ ip H (trx + tpH 



(26) 



Except for the equation V 4 tr% and V 3 tr%, the structure of the nonlinear terms in the other 
equations follow from signature considerations Q We now write the constraint equations (22) 



7 Notice that we have written a more precise version of schematic equation for V377 compared to V477. 
This will be useful in the proof since when integrating in the u direction using the V3 equation, we will not 
have a smallness in the length scale and we need to use the extra structure of the equation. 
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in schematic form: 

div X = ^Vtrx + ^(tr* + x)-f), 



div x = o Vtr X + ^( tr X + X) + §_ (27) 
curl 77 = —curl t] = a 
K = —p + ipip 



We now write down the Bianchi equations (24) in schematic form, substituting the Codazzi 
equations in (27) for some (3 and (3. In these equations, the left hand side is written with 
exact constants while the right hand side is written only schematically. 

V 3 /3 - Vp -* V<r =<Kp, <x) + ^ i2 (ipH + tr X )V i3 (V>H + trx), 

V 4 a + div */3 =^cr + ^ ^ V^V^Vff + ^XX, 

V 4 p + div/3=# + ^ ^V^V^h + ^XX, 

*l+«2+i3<l 

v 3 ct + div */3 + V' il v i2 V'V i3 VH + V'M' 

V 3 p + div /3 =# + ^ ^ V i2 ^ V i3 ^H + ^XX, 
n+«2+«3<i 

V 4 /3 + Vp -* V<r =^(p, a) + iJ H V h (^ H + tr X )V i3 (^H + trx), 

It is important in the sequel that in the equations for V 4 (p, a) (resp. Vs(p, a)), ipu (resp. 
ipn) does not appear. This does not follow from signature considerations alone since in 
principle the conservation of signature would allow a term i^h^h^h (resp. ^h^hJ^h)- The 
fact that these terms do not appear can be observed directly in the equation (24). 

2.6. Integration. Let U be a coordinate patch on S^o and pu be a partition of unity in Djj 
such that pu is supported in Du- Given a function 0, the integration on S UA is given by the 
formula: 



p poo poo 

/ := / / <&W det ldd 1 dd 7 

J S,,„ TT J —OO J —OO 



Su,u u oo J — CO 

Let -D u ' jtt ' by the region < u < u', < u < yf. The integration on D uu is given by the 
formula 

p pu pu ^00 poo 

/ : =5Z/ / / / #f/ V 7 - det gd6 1 d6 2 dudu 

J Du,u u J */ J — 00 J — oo 

=2 J2 / 0R/ft 2 V 7 - det ~fd9 1 d6 2 dudu. 

jj Jo JO J —00 J —oo 

Since there are no canonical volume forms on i7 u and H u , we define integration by 

P PE POO POO 

0-=J2 / 02p c/ fi v / det7^ 1 ^ 2 d«, 

1/ iJ u rr 1/ J —00 J —oo 
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and 

p re pco roc 

<t>-=J2 / <p2pu^&et-id6 l de 2 du. 

J H u jj Jo J— oo J — oo 

With these notions of integration, we can define the norms that we will use. Let <fi be an 
arbitrary tensorfield. For 1 < p < oo, define 

p 

V 
V 

L"(KJ 



:= f <<P 




J S u ,u 


:= f < 0, 




J H u 


= [ <0, 







Define also the L°° norm by 



\\<f>\\L~(s u , u ) := sup <0,0>f (5). 

. it 



We will also use mixed norms defined by 



L2L-LP(5) = ( / (SUP ||0|Up(S„,„)) <*U 

o ue[o,u»] 



LlL^LP(S) = / (SUP ||V>|| L p (5ii ,„)) Gfc/ 

o «e[o,e] 

Note that L°°L P is taken before taking L 2 . In the sequel, we will frequently use 

II ■ Wl^lzlp(s) < II • ||r»£«Li>(S). 

With the above definition, 1 10| \l 2 l 2 (s U!U ) an d | \4>\ \l 2 (h ) differ by a factor of f2. Nevertheless, 
in view of Proposition [TJ these norms are equivalent up to a factor of 2. 



2.7. Norms. We now define the norms that we will work with. Let 

'U 2 (SJ 



n = J2\ J2 su pII v ^IUw+ Yl su pIIv^I 

i<2 We{/3,/3,a} U *e{/5,a,/3} - 



ft(S) = ^(supHV^Aa,^)!!^^) + HV^II^^a^), 

. ... u.u 
1<1 '— 

= sup||V l (trx,r/,r/,trx)||Lf(5 u ,„) + ||V l (x, w)\\r%L°°i*(S) + l|V*(x, ^)Wlil^lp(s), 

u,u 

03,2 =||V 3 (trx,trx),V 2 (/i,/i, k,k)\\l^l^l 2 {S) + ||V 3 (^,^)|| LS o L 2 L2(5) + ||V 3 (^,?7)|| LflL 2 L 2 (s) 

+ ||V 3 (x,w,w t )|| L ^ L 2 i2(5) + \\V 3 (xM,^)\\l^i^(S), 
where and u} are defined to be the solutions to 

^ t 1 „ t 1 

V3W 1 = -a, V4W 1 = -cr 
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with zero data and fjj^ fJj^ K £1X6 defined by 
H := — div rj — p, ji := — div 77 — p, k := Vu +* Vw^ — -j3, K := -Vw +* Vu} — -(3. 



Recall from (25) that we use the schematic notation ip e {tr%, 77, r),trx}, e {x,^} 
and G {x, wj-T The choice of this notation is due to the fact that they obey different 
estimates. 

For the norms of the third derivatives of the Ricci coefficients, i.e., the Os^ norms, notice 
that V 3 tr% and V 3 tr% obey the same type of estimates as for lower order derivatives. V 3 (?7, rj) 
can no longer be controlled on a 2-sphere, but it obeys estimates on either null hypersurface. 
V^ipH (resp. V 3 ipn) satisfies similar estimates as before, but at this level of derivatives, we 
have to take L\ (resp. L 2 ) before (resp. L™). 
We write 

i<l i<2 

3. Statement of Main Theorem 

With the notations introduced in the previous section, we formulate a more precise version 
of Theorem |2j which we call Theorem [3] As noted before, since the proof for the existence 
and uniqueness of solutions in {0 < u < e} n {0 < u < I{\ is the same as that in {0 < u < 
e } H {0 < u < I 2 }, we will focus on the latter case. 

Theorem 3. Suppose the initial data set for the characteristic initial value problem is given 
on H for < u < and on Hq for < u < u* with < e and u* < I such that 

c < | det7 \ Su>0 |, | det7 \ So ^ \ < C, 

kol + l(^) 4 7 k„ \)<c, 

i<3 ^ ' 
°0 (ii V Vlk^(S u , ) + ||VVlU r ^(5o l34 ) + IIVVfflb(Ho) + l|V%ll^(Ho)) 

<c, 

i<2 y *e{p,d-} 
<C. 

Then for e sufficiently small depending only on C , c and I, there exists a spacetime (A4,g) 
endowed with a double null foliation u, u that solves the characteristic initial value problem 
to the vacuum Einstein equations in < u < I, < u< e. The metric is continuous in the 
(u,u, 6 ,1 ,# 2 ) coordinate system and takes the form In addition, 

g = -29? (du ®du + du® du) + ^AB{d6 A - b A du) <g> (d9 B - b B du). 

(A4,g) is a C° limit of smooth solutions to the vacuum Einstein equations and is the unique 
spacetime solving the characteristic initial value problem among all C° limits of smooth 
solutions. Moreover, 
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d 2 



86 

d 2 



2 geC° u ClL 2 (S), 



d6du 



g e L< U L™L\S) 



" 2 -g 6 LlLlL\S). 



dudu 

In the (u,u, O 1 ^ 2 ) coordinates, the Einstein equations are satisfied in L 2 U L 2 U L 2 {S\ Further- 
more, the higher angular differentiability in the data results in higher angular differentiability. 

In the remainder of this paper, we will prove the a priori estimates needed to establish 
Theorem [3] (see Theorem aprioriestimates) . The existence, uniqueness and regularity state- 
ments in Theorem [3] follow from the a priori estimates and an approximation argument as 
in [25]. We refer the readers to [25] for details. In the subsequent sections, we will prove the 
following theorem on the a priori estimates: 

Theorem 4. Suppose the initial data set for the characteristic initial value problem is given 
on H for < u < and on II for < u < with < e and w* < / such that 

c<|det 7 r Su , |<C7, J2\(^ls u , \<C, 

i<3 

°0 (H V >lk^ 2 (S u ,o) + IIVVIU^So,,) + IIVVffll^(Ho) + IIVVfdl^Ho)) 

<c, 

i<2 y *£{p,a-} 

<c. 

Then, there exists e depending only on C , c and I such that in < u < w*, < u < u# the 
following norms are bounded above by a constant C depending only on C , c and I: 

0,63,2,11 <c. 

3.1. Structure of the Proof. We briefly outline the proof of Theorem [4[ 

STEP 0: Assuming that C\),oo and Ci i4 are controlled, we prove the bounds on the metric 
components, from which we derive preliminary estimates such as Sobolev embedding and 
the estimates for transport equations. (Section El). 
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STEP 1: Assuming K < oo, TZ(S) < oo and (9 3i2 < oo, we prove that O < C(O ,K(S)). 



(Sections 5.1 , 5.2) 



STEP 2: Assuming K < oo and C 3 , 2 < oo, we show that K(S) < C(TZ ). (Section |5_3j) 
Together with Step 1 this implies O < C{O ,n ). 

STEP 3: Assuming K < oo, we establish that 3 , 2 < C(O )(l + K), i.e., 3 , 2 grows at 



most linearly with 1Z, with a constant depending only on the initial data. (Section 5.4) 



STEP 4: Using the previous steps, we obtain the estimate TZ < C(Oo,TZo), thus finishing 
the proof of Theorem 111 (Section [6j 



4. The Preliminary Estimates 
All estimates in this section will be proved under the following bootstrap assumption: 



0O,oo + $^i,4< A . (Al) 



t<l 



4.1. Estimates for Metric Components. We first show that we can control Q under the 



bootstrap assumption (Al): 



Proposition 1. There exists eo = e (A ) such that for every e < e , 

1 



< n < 2. 

2 ~ ~ 



Proof. Consider the equation 



1 1 X c) 

oo = --V 4 logfi = -IWT 1 = 2^ irl - ( 28 ) 

Fix u. Notice that both u and Q are scalars and therefore the L°° norm is independent of 
the metric. We can integrate equation (28) using the fact that Q^ 1 = 1 on if to obtain 

f~ i 

_ l|U°°(Su,«) ^ C I \\ UJ \\L°°(s uu ,)dyf < Cez 1 1^11^1,2^(5) < CA e2. 
Jo 

This implies both the upper and lower bounds for Q for sufficiently small e. □ 



We then show that we can control 7 under the bootstrap assumption (Al): 



Proposition 2. Consider a coordinate patch U on Sq } o and define U Ui o to be a coordinate 
patch on S U) o given by the one-parameter diffeomorphism generated by L. Define U u ^u to 
be the image of U Ui q under the one-parameter diffeomorphism generated by L. Define also 
Du = Uo<n</ o<«<€ ^u,u- For e small enough depending on initial data and A 0; there exists 
C and c depending only on initial data such that the following pointwise bounds for 7 in V>u 
hold: 

c < det 7 < C. 

Moreover, in Djj, 

\1ab\, \(l~ 1 ) AB \ < C. 
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Proof. The first variation formula states that 

In coordinates, this means 

d 

— 7ab = Z^Xab- 
ou 

From this we derive that 

d 

— log(det 7) = fttrx- 
Define 7o (u, u, 1 , 2 ) = 7(1*, 0, 1 , 2 ). Then 

I det7 - det(7o)| < C / |trx|du' < CA e. (29) 

Jo 

This implies that the det7 is bounded above and below. Let A be the larger eigenvalue of 
7. Clearly, 

A<C7 sup 7ab, (30) 

A,B=1,2 

and 

\Xab\ 2 < CA\\x\\l°°{s u ,u)- 

A,B=1,2 

Then 

\1ab ~ {iq)ab\ <C \xab\oIu < C7AA e5. 
^0 

Using the upper bound (30), we thus obtain the upper bound for |tab|- The upper bound 
for |(7~ 1 ) AB | follows from the upper bound for |7ab| and the lower bound for det7. □ 

A consequence of the previous proposition is an estimate on the surface area of the two 
sphere S u& . 

Proposition 3. 

sup \Area(S U: u) — Area(S u<0 )\ < CA e. 



Proof. This follows from (29). □ 

With the estimate on the volume form, we can now show that the LP norms defined with 
respect to the metric and the L p norms defined with respect to the coordinate system are 
equivalent. 

Proposition 4. Given a covariant tensor (pAi...A T on S UtU , we have 

! <0,0>f~]T Y, If ^...A^V^fde'de 2 . 

JSu.n i=i Ai=l,2 

We can also bound b under the bootstrap assumption, thus controlling the full spacetime 
metric: 

Proposition 5. In the coordinates system (u,u, 6 l ,6 2 ), 

\b A \ < CA e. 
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Proof. b A satisfies the equation 



du 

This can be derived from 

[LM d " A 9 



f)h A 

-4tf( A . (31) 



du d6 A ' 



Now, integrating (31) and using Proposition 4 gives the result. □ 



4.2. Estimates for Transport Equations. The estimates for the Ricci coefficients and the 
null curvature components are derived from the null structure equations and the null Bianchi 
equations respectively. In order to use the equations, we need a way to obtain estimates 
from the covariant null transport equations. Such estimates require the boundedness of 



tr% and tr%, which is consistent with our bootstrap assumption (Al). Below, we state two 
Propositions which provide LP estimates for general quantities satisfying transport equations 
either in the e 3 or e 4 direction. 

Proposition 6. There exists €q = eo(Ao) such that for all e < €q and for every 2 < p < oo, 
we have 

pu 

LP(S u ,u) <C(\\(j)\\LP(S uul )+ / | |V 4 0|Up(5 u „/,)<*!*")> 

Ju> 
r-u 

LP(S Ul3i ) < C(||0||if(5 tt , J + / ||V 3 0|| iP( 5 u „ u )du"). 

Ju> 

Proof. The following identity holds for any scalar /: 



Similarly, we have 

d 



I f = t n (e 3 (f) + trxf) . 

S11..11. J St 1, .11, 



du 

Hence, taking / = \<j)\Z, we have 

1 M \Us^j =M I Us uvl ) + [~ f Pl^r 2 ^ (< 4>> V 4 > 7 +^tr X |0|2 ) du" 



P 



[ U [ pW~ 2 n (< 0, V 3 > 7 +-trx|0|^ du'' 

Ju> Js„„„ \ P J 



(32) 



The proposition can be concluded using Cauchy-Schwarz on the sphere and the L°° bounds 
for n and tr% (tr%) which are provided by Proposition [I] and the bootstrap assumption (Al) 
respectively. U 

The above estimates also hold for p = 00: 

Proposition 7. There exists €q = e (A ) such that for all e < 6q, we have 

|l°°(s„, h ) < C ( ||0lU°°(s u „,) + / ||V 4 0|U~(s uu „)tft/ 



{Su>%) < C \\4>\\l°°{s u , u ) + / ||V 3 0|U~ (5 jdu" . 
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Proof. This follows simply from integrating along the integral curves of L and L, and the 
estimate on Q in Proposition [Tj □ 

4.3. Sobolev Embedding. Using the upper and lower bounds of the volume form, Sobolev 
embedding theorems in our setting follows from standard Sobolev embedding theorems (see 

!):• 



Proposition 8. There exists €q = eo(A ) such that as long as e < €q, we have 



i=0 

Similarly, we can also prove the Sobolev embedding theorem for the L°° norm: 
Proposition 9. There exists e = e (A ) such that as long as e < e , we have 

UWl^s^) < c (UWl^) + HV0|| L 4 ( ^ m) ) . 

As a consequence, 

2 

Besides the Sobolev embedding theorem on the 2-spheres, we also have a co-dimensional 
1 trace estimate that controls the L 3 (S) norm by the L 2 (H) norm with a small constant. 

Proposition 10. 

IMIls(s WiH ) <C (lHU 3 (s u ,„0 + ^||V0|| L 2 i2(5) + ei||V 4 0||L2 L 2 (5) ) . 

Proof. It follows from standard Sobolev embedding theorem and the lower and upper bounds 
of the volume form that 

< c(!|0|i! 3(5) ||V0|i! 2(5) + M\ LH s)). (33) 



Using (32) and (33), we have 

M\h iSu , u) =M\Us uul ) + T f 3fi|0| 7 (< 0, V 4 > 7 +^tr X |0|2 ) dyf' 

Ju' J s„, V 6 



<M\h(,s u u ,) + m\l H H) 1 1 V 4 0| \l*( H ) + f CA 1 10| \1 3{S jdu" 

Jo 

+ rcA M\h iSuu/) dyf 

Jo 

<M\h^) + (M\l~ LH s)\\Vm kLHS ) + ^ll0llig=L3 ( 5))l|V40|U 2 (H) 

+ rC7Ao|H|i, (5 jdu'. 
Jo 
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Using Holder inequality and absorbing the term 1 10| ||oo L 3f S ) to the left hand side, we have 



3 



<c [M\h {Su , u _,) + IIV0IU^(5)l|V40||i| L2(5 ) + ^llv 4 0|li iL 2 {5) 

+ CAoMUs^ uf) dy! 
Jo 

<C (ll^llia^,) + ^||V0|| 3 LiL2(5) + el||V 4 0||i iL2{5) +ei||V 4 0||l liS 



L3(S) 



By Gronwall's inequality, and using the fact that u < e, we have 

UWl^) <^(ll0lli3(5^)+ el ll V ^li|L 2 (5)+ el ll V 40lli|L 2 (5))- 

□ 

4.4. Commutation Formulae. We have the following formula from [18] : 
Proposition 11. The commutator [V 4 , V] acting on an (0,r) S-tensor is given by 
[V 4 , V B )(f> Al ... Ar =[D 4 ,D B )(f> Al .., Ar + (V B logn)V 4 0A 1 ...A r - (T 1 ) CD XBD^c<pA,...A r 

r r 

~ J2^^ CDxBDr lA^A 1 ...A l c...A r + J2^~ 1 ) CD XA i Br[ D <f> Al ^ A .a„ Ar . 

i=l i=l 

By induction, we get the following schematic formula for repeated commutations (see [25J): 
Proposition 12. Suppose V 4 = F . Let V 4 V 4 = F,. Then 

il+«2+«3=« h+i2+i3+ii=i 

where by V H (r] + rf) 12 we mean the sum of all terms which is a product of i 2 factors, each 
factor being V J (r] + r() for some j and that the sum of all j's is i\, i.e., V ll (r] + rf) 12 = 

V J1 (?7 + 7])... V Ji2 (rj + ri). Similarly, suppose = Gq. Let VsV*0 = G{. Then 

ji+-+ji 2 =h 

Gi ~ Yl v h (v+y) i2 v i3 G + v h (v+vy 2 v i3 xy i4 

U+«2+*3=« jl+«2+«3+«4=i 

+ V il (r] + r L ) i2 V i3 §V i4 (f). 

«l+«2+«3+*4=«— 1 

The following further simplified version is useful for our estimates in the next section: 
Proposition 13. Suppose V 4 = Fq. Let V 4 V*0 = Fi. Then 

Fi~ VV 2 V 43 F + V^^V^xV* 4 ^ 

il+h+i3=i U+i2+«3+*4=i 
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Similarly, suppose V30 = Go. Let VaV l = Gi. Then 

*l+«2+«3=* il+i2+*3+*4=« 

Proof. We replace /3 and /3 using the Codazzi equations, which schematically looks like 

P = Vx + i/)X, 
P = Vx + *Px- 

□ 

4.5. General Elliptic Estimates for Hodge Systems. We recall the definition of the 
divergence and curl of a symmetric covariant tensor of an arbitrary rank: 

(div (j))A 1 ...A r = V B (j)BA 1 ...A r , 

(curl 4>) Al ...A r = { BC V B (ficA 1 ...A r , 
where j. is the volume form associated to the metric 7. Recall also that the trace is defined 
to be 

(tr(f>) Ai...A r -i = {T l ) BC <i>BCAr...Ar-x- 

The following elliptic estimate is standard (See for example [H] or [8]): 

Proposition 14. Let be a totally symmetric r + 1 covariant tensorfield on a 2-sphere 
(S> 2 ,7) satisfying 

div 4> = f, curl (ft = g, trxj) = h. 

Suppose also that 

^||V^|| L 2(5) < OO. 

i<l 

Then for i < 3, 

i-l 

IIVVIU^S) < C(%2 HV fc ^|| L2 (5))E(||VVI|L 2 (5) + \\V j 9\\LHS) + l|V^|| L2( 5) + MIl^S)))- 
k<l j=0 

For the special case that a symmetric traceless 2-tensor, we only need to know its 
divergence: 

Proposition 15. Suppose <fi is a symmetric traceless 2-tensor satisfying 

div 4> = f. 

Suppose moreover that 

^||V^|| L 2(5) < OO. 
i<l 

Then, for i < 3, 

i-l 

||VV|U" ( 5) < C(£ \\V k K\\ LHs) )(Y^(\\^f\\ LHS ) + MIlhs)))- 

k<l j=0 



Proof. In view of Proposition 14, this Proposition follows from 

curl (ft =* f. 

This is a direct computation using that fact that is both symmetric and traceless. □ 
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5. Estimates for the Ricci Coefficients 



We continue to work under the bootstrap assumptions ( Al ). In this section, we show that 
assuming the curvature norm 1Z is bounded, then so are the Ricci coefficient norms O, 0^,2 



and the curvature norm 1Z(S) on the spheres. In particular, our bootstrap assumption (Al) 
and all the estimates in the last section are verified as long as TZ is controlled. 



5.1. L 4 (5) Estimates for First Derivatives of Ricci Coefficients. 
Proposition 16. Assume 

1Z < oo, (5 3)2 < oo, 2 ,2 < oo. 
Then there exists e = e (7Z, 3)2 , 02,2, A ) such that whenever e < e , 

Y,OiAtrX,v]<C(Oo). 

i<l 

In particular, C(Oq) is independent of A . 

Proof. Using the null structure equations, we have a schematic equation of the type 

V 4 (trx, v) = ft + p + Vn + ipij) + ipHip- 

It is important to note that (3, if)jj do not appear in the source terms. In other words, only 
the terms that can be controlled on the outgoing hypersurface H u enter the equation. We 
also use the null structure equations commuted with angular derivatives: 



»l+i2+»3=* »l+»2+»3-H4=» 

By Proposition 6, in order to estimate ||V l (tr\;, \l°°l°°l 4 (S), it suffices to estimate the 
initial data and the || • \\l°°l^l a {S) norm of the right hand side. We now estimate each of the 
terms in the equations. For the curvature terms, we have 

|| ^ i2 (^P)\\L-L k LHS) 
ii<l «2<1 

<C(l + A )6^||V i (/3,p)|| Lsoi|L2{5 ) 

i<2 

<C(1 + A )e^K. 
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The term with V77 instead of (/3, p) can be bounded analogously, except for using the and 
norms together instead of the TZ norm: 



|| J2 ^V^vWl-lilhs) 

ii<l i 2 <2 

<C7(1 + A )e5 ^||V i r/|| L ^ LrL 4 (s) + ^ ||VS,||^ ( J 

\i<l 2<i<3 / 

<C(1 + A )e3 O iA + 2 , 2 + 3>2 ) 



»<i 

1 



<L7(l + A )e2(A + 2i2 + 3 , 2 ). 
We now move on to the lower order terms: 

11 22 ^vh^hwl^lhs) 
<c(E iMiW- ( 5))(E ii vi ^iu^(s)) 

il<2 i 2 <l 

<L7A (1 + A ) 2 e. 
Finally, we bound the lower order terms that contain ■?/>#: 

|| ^ ^V^W^IU-^s) 

<c(E IMi W°W(E liv^ll^^w) 

il<2 i 2 <l 

+^(Ei^iiW^w)(Eii via ^ii^ £ ^))(ii^iu-i i x«(5); 

ii<l i 2 <l 

<C(1 + A ) 2 e5(^ ||Vtyir|U»La^(s) + E ll V VIUw^(S)) 



i<l i<l 
1 



<0A Q (1+A O )^ 
Hence, by Proposition |6j we have 

22 Oi,i[trx, V] <&o + C(l + A ) 2 e3 (K + 2 , 2 + 3 , 2 + A ). 

The Proposition follows from choosing e to be sufficiently small, depending on 71, 03,2, 02,2, A . 

□ 

We now estimate the terms that we denote by ipn, i-e., x an d Both of them obey a V4 
equation. However, a new difficulty compared Proposition 16 arises since the initial data for 
X and uj_ are not in L£°. Thus they can only be estimated after taking the l? u norm. 

Proposition 17. Assume 

1Z < OO, 03 2 < OO, 2 2 < OO. 
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Then there exists e = €q(TZ, 3>2 , 2j2 , A ) such that whenever e < e , 



i<l 



In particular, this estimate is independent of A . 

Proof. Using the null structure equations, for each ipH £ we have an equation of 

type 

W 4 ipH = P + Vr/ + (ip + ^h){^ + ^h)- 
We also use the null structure equations commuted with angular derivatives: 

V 4 VV^= Yl VV 2 V 43 (p + Vr?) + J2 V^V^ + ^V^ + V^) 



From the proof of Proposition 16 , we have 



V Jl V i2 V i3 p|| Ls ,^L4 (5) <C(1 + A )e^, 

*l+*2+«3<l 



and 



^ V i2+1 r7| 1^^4(5) <C(1 + A )e5 (A + (9 2 , 2 + d 3 , 2 



«1+*2<1 



and 



£ ^^VV^VlUyx^^ <CA (1 + A ) 2 e, 

«l+*2+«3<l 



and 



^V^V^hWl-lilhs) <CA (1 + A ) 2 e3 



U+«2+«3<l 



The two new terms that did not appear in the proof of Proposition 16 are 



V* 1 V> <2 V is VirV <4 Vir, and £ V n # 2 V^V^a- 

U+«2+«3+*4<l jl+«2+«3+M<l 
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Both of these terms cannot be controlled in the L^L^L^IS) norm. Instead, for each fixed 
u, we bound the first term in the L^L 4 (S') norm: 

n+«2+«3+*4<l 

<C(1 + \\iP\\l^L^(S))(\\^h\\l^L^(S))(Y^ I |VVtf| I L^£4(S)) 

i<l 

+ C(l + \\lp\\ L ™L°°(S))Q2 II V * 3 V ; hIU^L 4 (5))(||V ; h|U^L-(S)) 

i<l 

<C(1 + A ) 2 e5(||^||L2L-(s) + ^\\V1>h\\izlHS)){^\\^h\\l?l*(S)) 

i<l i<\ 

<c(i + a )M(^ ||vVdU|^(5) + X) II v VhIU^ 4 (5))- 

i<l i<l 

by Sobolev Embedding in Proposition [8] 
<C(1 + A )M(A + ^ HVVslU^^). 

i<l 

According to the definition of the norm, ^ obeys stronger estimates than Therefore, 
we can control the remaining term in the same manner: 

1 1 J2 VH ^ 2 ^H^Hh\ \ L ims) < C(l + A ) 2 e5 (A Q + ^ || VtyflJ |^(S)). 

Therefore, by Proposition [6j for all u G [0, e], 

J3ll V ^IU 4 (s«,«) 

i<l 

<C£ HVValU*(S^o) + (! + Ao) 2 e5(^ + 2 , 2 + 03,2 + A + IIWh|U»x*(S)))- 

i<l Kl 

Clearly the right hand side is independent of w. Thus we can take supremum in u on the 
left hand side. Then, we take the L\ norm to obtain 

^||VVfdlz^£4(S) 
i<l 

<C(^ ||VVgJ|^(s Ul0 ) + (1 + A ) 2 ^ (ft + 3 , 2 + A + ^ 1 1 VVjd 

i<l i<l 

<c7(O + (1 + A ) 2 e5(^ + 2 , 2 + 3 , 2 + A )). 
The left hand side is precisely what we need to control for the Ci,4[X)^] norm. Thus 

i<l 

°iA%_M < C(&o + (1 + A )M(ft + 2 , 2 + a 3 , 2 + A )). 

i<l 

We conclude the proof by choosing e to be sufficiently small. □ 

We now turn to the Ricci coefficients T], x, u. To estimate these Ricci coefficients, we use 
the V 3 equations. Unlike in the proofs of Propositions [16] and [TT| where a smallness constant 
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can be gained from the shortness of the u interval, when integrating the V3 equation, the 
u interval is arbitrarily long. Instead, we show that the inhomogeneous terms are at worst 
linear in the unknown and the desired bounds can be obtained via Gronwall's inequality. 
Notice that x satisfies a V4 equation with a as a source term. We avoid this equation 
because a is singular. We begin with the estimates for 77. As we will see below, we cannot 
directly estimate its L 4 (S) norms, but have to first estimate its L°°(S) norm: 

Proposition 18. Assume 

1Z < OO, O3 2 < OO, 02,2 < 00. 

Then there exists e = €q(TZ, (^3,2, ^2,2, A ) such that whenever e < e$, 

Oo,oo[v]<c{o o ,n{s)[0\). 

In particular, this estimate is independent of Aq. 

Proof. 77 satisfies a V3 equation. Integrating in the u direction does not give a small constant 
as in integrating in the u direction. We therefore need to exploit the structure of the equation. 
We have, schematically 

V 3 ?7 = (ipH + tr%) (77 + 77) + /?. 

We notice that the quadratic term rf does not appear. Moreover, tr%, \ and u> do not 
enter the equation. In other words, all Ricci coefficients except 77 in this equation have been 
estimated in the previous propositions by C(O ). We now bound each of the terms. Firstly, 
the term with curvature can be controlled using Holder's inequality and Sobolev embedding 
in Proposition [9] by 1Z(S): 

\\P\\l~lil~w < Cj^W^tW^Lims) < C^R(S)[§]. 
Then, we estimate the terms quadratic in the Ricci coefficients, which do not involve 77: 

\\{lpH + tlx)v\\L^LlL^(S) 

<C\\ipH + ^x\\l^lil^(s)\\v\\l^l^l^(s) 
<C{O ), 

by Propositions 16 and [17] and Sobolev Embedding in Proposition [9} Finally, we estimate 
the term (ipu + trx)?]- Fix u. Then 

\\{ipH + trx)v\\ L i L ™ iS ) 

P U 

<C / \\if)H + trx\\Lo°(s u , J\v\\loo { s u , jdu'. 
Jo 

Therefore, by Proposition [6j we have, for every u, 

\\v\\L?i^(S)<C(Oo) + CI*R(S)\g+C / ||^ + trx|U- ( s u ,j|N|L»(s u , jdu'. 

Jo 

By Gronwall's inequality, we have, for every u, 

\\v\\l~l°°(s) < C{O ,n{S)[§])e W {C [ \WH + tTx\\L~(s n ,jdu'). 

Jo 
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Using Cauchy-Schwarz, the Sobolev embedding in Proposition |9l as well as the estimates for 



the Ricci coefficients derived in Propositions 16 and [TTJ we have 

\\v\\l~l~(S)<C{O q ,TI{S)\1]), 
as desired. □ 

Using the L°° estimate of 77, we now control V77 in L 2 : 
Proposition 19. Assume 

1Z < 00, d) 3)2 < 00, 2 ,2 < 00. 
Then there exists e = €q(TZ, (^3,2, ^2,2, A ) such that whenever e < e$, 

O 1>2 [r]}<C(O o ,n(S)[0\). 
In particular, this estimate is independent of Aq. 
Proof. Recall that we have, schematically, 

V35 = {^h + trx) (77 + 77) + p. 
Commuting with angular derivatives, we get 

v 3 V7/= J2 (v + y) h v i2 f3_+ J2 (v + y) h V i2 (v + v)V i3 (iJH + trx)- (34) 

il+i 2 =l ii+i2+«3=l 



We notice that in (34), when two 77's appear in a term, only one of them has a derivative. 
Fix u. We now estimate each of the terms. Firstly, the term with curvature: 

|| Yl (v + vY^PWl^s) 

<c(i + 11(17, ^^^(^(x; ii v ^iu^^ 



i<l 

i<l 

<C(O ,K(S)[P\), 



since V/3 can be controlled in L 2 (H_ U ) by TZ(S). We then estimate the nonlinear term in the 
Ricci coefficients: 

|| (v + vy* V i2 (77 + 77) V 13 (Vh + trx) 1 1^(5) 

il+«2+i3<l 



< 



C(l + 1 1 (77, 77) ||l-l-(5)) 2 ($^ 1 1 V 4 (V s , trx) I kL»(fl)) 



i<l 



+ C7 / ||V(7y,5)||£3 (Su; j 1 1(^, trx) 
<L7(O , [£]) + C [ U 1 1 V77J \ LHSu ,, j 1 1 trx) ||L»(s„, ifi )du', 
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where the first term is bounded using Propositions 16 17 and 18 Therefore, by Proposition 
[61 we have, for every u, 



V\\lz°l 2 (Su,u) 



<C(O ,n(SM) + C / \\V R \\ LHSu ,M^K^x)\\L-(s u , u )du'. 

Jo 



By Gronwall's inequality, we have 



<C(O 0} K(Sm)eM / \\(i>H,trx)\\L°o(s u , jdu'). 

Jo 



The right hand side satisfies the desired bound by Propositions 16 and 17 



□ 



Recall that by Proposition 18 we now have a bound on Oo,oo[^] independent of Aq. This 



allows us to prove the Oi^^jestimates. However, unlike the L°°(S) estimates for r\ and 
L 2 (S) estimates for V//, the L A (S) control that we prove at this point for V77 grows linearly 
in TZ. This bound will be improved in the next subsection. 

Proposition 20. Assume 

TZ < OO, 03,2 < OO, 02,2 < OO. 

Then there exists eo = eo(7Z, ^3,2, ^2,2, A ) such that whenever e < eo, 

O 1A [n\<C(O ,K(S))(l + K). 
This estimate is linear in the TZ norm and is independent of Aq. 
Proof. Recall that we have, 

V 3 V^= J2 (v + v) h V i2 f3_+ J2 (v + v) h V i2 (ri + ri)V i3 (tlJH + trx). 



ii+i2+»3=l 



As in the proof of Proposition 19, we notice that in this equation, when two 77's appear in a 
term, only one of them has a derivative. Fix u. Now, we estimate each of the terms. Firstly, 
the term with curvature: 



< 



c(i+\\(v,m^ {s) )(j2\\^t\\^w 



i<l 



< 



c(0 o ,rc(s))X)ll v< £lk*» 



(S) 



i<2 



<C{O ,K{S))K. 
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We then control the nonlinear term in the Ricci coefficients: 

|| J2 ( V + 77)^(77 + 77)V i3 (^ + trx) I \ K lHS) 



< 



ca + ii^^iu^^^^iiv^^trx)!!^ 



i<l 



+ C / \\V(r),ri)\\L H s u , J\(i)H^*x)\\L°°(s u , jdu' 
Jo 

<C(O , K(S)) + C T \\Vrj\ \ LHSu , j 1 1 (fe tax) I U~(s u , „)**'. 



Therefore, by Proposition [6j we have, for every m, 

\\Vl±\\ L ce LH S) 

<c(o ,n(S))(i + n) + c P \\Vv\\ LHSu , Jfetrx)!^, j*/. 

^0 

By Gronwall's inequality, we have 
||Vr7|| L ^ L 4 (5) 

<C(Oo,^))(l+^)exp(r||(fetrx)||Lo 0(Su , J ^)- 

./o 

The right hand side satisfies the desired bound by Propositions 16 and 17 
We now estimate the norm of ipn- 

Proposition 21. Assume 

1Z < OO, 3j2 < OO, 2 ,2 < OO. 

TTien there exists e = e (7?., 03,2, ^2,2, A ) s,uc ^ whenever e < e , 

J]Om[x.w] <C(O O ,K(S)[0\). 

i<\ 

In particular, this estimate is independent of A Q . 
Proof. Consider the following equations for ipn £ {x,u}: 

^3^H = V?7 + p + ip H (trx + ^h) + 4>(4> + i>H)- 
As before, we commute the equations with angular derivatives: 



h +«2=1 



+ (^ + !Z) il (V i2 ^V i3 (trx + ^) + V i2 ^V i3 (^ + ^//))- 



U+i2+i3+*4 = l 
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We bound each of the terms in L l u L 4 (S). First, we look at the curvature term: 



<C(1 + \\(r),v)\\Lz>L°°(S))(Y^ \\V 1 p\\l~lil±(s)) 

i<l 

<C(A )J2\\V*p\\ LlL 2 {s) 

i<2 

<c(A ,n). 

The term containing V 2 ry can be estimated analogously: 



ii (v+v) h ^ +x v \\ LiLHS) 

<c(i + \\(v,v)\\ L ^(s))(Yl Hv i+1 ?7lU r iiL4 (s) ; 

i<l 

<C(A )J2\\V i+1 v\\LlL>( S ) 

i<2 

<c(a q ,omo 2 ,M)- 



Then, we control the terms containing ipn'- 



j2 (v+vy^^H^^H+tix)^ 



(S) 



pu 

<C\\(v,v)\\l^l^(s) / \\iPh\\l^(s u , j^HVX^trx)!! 

JO ' ^ 



du 



i<l 



+ C\\(v>v)\\l°°l<*>(S) 



< 



/ J2 1 1 v Vir| U*(s„,,„) 1 1 tT X) \\L°°(s u , !U )du' 

<7(0 o ,rc(£)[g]) r^iiv^^ll^^.j^liv^^trx)!!^^,,)^'. 

^° ii<l i 2 <l 



In the above, we noticed that and r\ obey estimates from Propositions 16 and 19 that 
depend only on Oq and 1Z(S)[(3}. For the terms not containing ipn, we can bound directly 



using the bootstrap assumption (Al) 



il<l «2<1 *3<1 

<C(Ao). 
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Therefore, by Proposition [6j 



i<l 



<C l|VVif||L4(5 , a ) + C(K, 6 3j2 [r)}, 2j2 [rj], A ) 

Jo ' • . , 



il<l «2<1 

By Gronwall's inequality, we have 



i<l 



<C(J2 \\V^h\\lhs 0a ) + C(K, 6 3 , 2 [ V ], O2M Ao)) 

x exp( I" C(O O ,K(S)[0\) J2\\V\^H,trx)\\ L Hs u ,jdu'). 



i<l 



By Propositions 16 and 17, we have 



exp( / L7(O ,^(S)[^])^||V J (^,trx)||L 4 (s u , fi) d«' < C{O ,K(S)\g). 



i<l 



Therefore, we have, for any u, u, 

5^ l|VVirlU*(Su^) 
»<i 

<c(0o,ft(£)[£])(X;^ 

i<l 

Clearly the right hand side is independent of u. We first take supremum in u and then take 
the L\ norm to obtain 

i<l 

1 

<C7(O , ^WI^DCX) HW^IIza^Cflb^ + c'C7(7i, <9s^W, 02,aW, Ao)). 

i=0 

By choosing e sufficiently small depending on 71, 3>2 , 2>2 , A , we have 

^||VVh||l|^^ ( 5) < c(0o,ft(S)[£]). 

i<l 

□ 

We then estimate the £^4 norm of tr%. Although tr% satisfies a V4 equation, the 

i<l 

term XX appears on the right hand side and each of the x factor has to be estimated in L\. 
Therefore, the bound for this term does not have a smallness constant. 
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Proposition 22. Assume 

1Z < OO, (5 32 < OG, C?2,2 < OO- 

T/ien i/iere exists e = e (72., (9 3) 2, 2t 2, A ) sitc/i £/iat whenever e < e , 

J2°iAtrx]<C(O o ,TZ(S)[0\). 

i<l 

In particular, this estimate is independent of Aq. 

Proof. Using the null structure equations, we have an equation of the type 

V 4 trx = tpip + 4>h^ + XX- 
We also have the null structure equations commuted with angular derivatives: 

V 4 V i trx= vV 2 v i3 ^v i4 (v> + <M + ^ h (V:V) i2 V i3 xy U X- 

ii+i2+i3+i4=i il+«2+«3+M=i 



By Proposition 
data and the 1 1 • 



6l in order to estimate IIVVIIl^z 00 !^^ it suffices to estimate the initial 
l™l\l a (s) norm of the right hand side. Notice that all terms except the one 



with XX have appeared in the Proposition 16 We estimate those terms in the same manner. 
Hence, 

|| ^V'W^IU-l^s) <CA (1 + A ) 2 e, 

and 

|| ^ i2 ^ i3 M\L~LlL*{ S ) <CA (1 + A ) 2 ei 



For the term with xXi using the estimates obtained in Propositions [19] and [2TJ we have 

U+i2+«3<l 

<C(J2 \\^^)\\l^l^l^( S ))(\\x\\l^lIl^(s))(J2 1 1^X11^^(5)) 

ii<2 «2<1 

<C(O ,K(S)[P}). 



Hence, 



J2°iA^x] <C(O O ,K(S)[0\) + CA (1 + A ) 2 ei 



i<l 



The proposition follows by choosing e to be sufficiently small depending on Aq. □ 



Clearly Propositions [16| [17j, |20| |21] [22j imply the following estimate for the L norms of 
the Ricci coefficients: 

Proposition 23. 

7Z < oo, 1Z(S) < oo, 3:2 < OO, 2 ,2 < oo. 
There exists e = e (C\), TZ, 7Z(S), 3 2 , ^2,2) such that for all e < e 0; we have 

i<l 
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J2OiAtrx,X,u;}<C{O o ,n(S)[0}), 



i<l 



and 



<C(O ,K(S),K). 



i<l 



Together with Sobolev embedding in Proposition^ the bootstrap assumptions (Al) can be 
improved under the assumptions on TZ,TZ(S), 03,2 and (Dis- 



proof. Let A > max{C(O ),C(O ,n(S) [/?]), C7(0 O , n{S),K)}, where the right hand side 



small so that the conclusions of Propositions 16 



171 [20 



21 


[22] 


21 


[22] 



embedding from Proposition [9j we have improved ( Al ). Since the choice of A depends only 
on Oo,lZ,7Z(S), the choice of eo depends only on O , 71, 7Z(S), 03,2, 02,2- □ 

5.2. L 2 (S) Estimates for Second Derivatives of Ricci Coefficients. We now estimate 
the 02,2 norm. We make the bootstrap assumption: 

2 ,2 < Al (A2) 
The proof of the estimates for the 02,2 norm is very similar to that for the 0j,4 norms, 

i<\ 

except that we now need to use the L 4 control that was obtained in the previous subsection. 
From now on, we will assume e < eo as in Proposition 23, where eo depends on Oq,TZ, 03,2 
and also on Ai. 

We first prove the estimates for V 2 tr% and V 2 ^: 

Proposition 24. Assume 

TZ<oo, 3j2 < 00, K(S) < 00. 

Then there exists eo = €o(Go,H, 03,2, 7l(S), Ai) such that whenever e < eo, 

2 ,2fer/]<C7(0 o ). 
In particular, this estimate is independent of A\. 
Proof. Using the null structure equations, we have 

V 4 (trx, rj) — P + P + Vr/ + vpt/j + tpHip- 
We use the null structure equations commuted with angular derivatives: 

V 4 V i (trx,r7) = J2 V n ^ 2 V is (/3 + p + Vr/) + ^ V n # 2 V i3 ^V i4 (^ + %j} H ). 



il+*2+«3=« 



*l+*2+«3+*4=« 



By Proposition 6, in order to estimate \ \V l {tTXi r i)\\L%>Lz>L 2 (s)i h suffices to bound the initial 
data and the || • \\l^l^l 2 (s) of the right hand side. We now estimate each of the terms in 
the equation. We first control the curvature term. As mentioned in the beginning of this 
subsection, the bounds are derived similarly as that for the L 4 norms, except we now need 
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to use the L 4 (S) estimates proved above for Vip. 
|| Yl V^ i2 V i3 (/3,p)|| i$fi i i2(5) 

<C(E IMlW"-(*))(£ \\VH(3,P)\\L~L k L HS) ) 
h<3 i 2 <2 

ii<2 i 2 <2 i 3 <l 

i<l 



i<l i<2 i<2 

l 



The term with Vi] instead of curvature can be estimated analogously, except for using the 
02,2 and 03,2 norms instead of the 1Z norm. Moreover, recall that the ^^£^4 [77] bounds 

Kl 

depend on 1Z. 

ii<3 j 2 <2 

+ C(E II^IU-L-L-(5))(E HV i2 ^||LWL2(5))(E H V13+1 ^IUW 4 (S)) 
ii<2 i 2 <2 i 3 <l 

<C(O ,^))6^((El|VVIU S =L r L^S)? + E ll V VIU-L^L 2 (S))(E H V ^HiS°iIi 2 (5)) 
Kl i<2 j<3 

<C(Oo,^(5),A 1 ,a 3 , 2 ,^)e5. 

We now move on to the lower order terms. We first control the lower order terms that 
contain ip H : 



Y V^V'W^IIl-l^s) 

n+*2+*3+*4<2 

<C^(E II^H W^(5))(E W^hWl-lil^s)) 

U<3 i 2 <2 
+ C^(E H Vn ^IU$f^i 4 (S))(E H Vi3 ^IU-L|L 4 (5)) 

U<1 1 2 <1 

+ Ce5(^ ||VlliVL^(S))(E l|V l2 ^||L r L 2 L 2 (5))(||^||L $f L|^(5)) 
U<1 i 2 <2 

<C(0 O , 7e)e5 (J2 1 1 II^l2l2(5) + 



<C(O ,7e(S'),7e,A 1 )e5. 



i<2 i<2 
1 
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The remaining lower order terms can be estimated in the same way since the norms for ip 
are stronger than those for ipu 

|| Yl VV 2 V*W^IIl-l^(s) <C(O ,n(S),TZ,A 1 )e 1 K 

«l+«2+«3+«4<2 

The conclusion thus follows from the above estimates and Proposition |6j after choosing e to 
be sufficiently small depending on O , 7Z(S), A 1; 3)2 , H- D 

We then estimate V 2 ipn_, where ipH_ G {x^}- 
Proposition 25. Assume 

TZ < oo, 3j2 < oo, K(S) < oo. 
Then there exists e = € (O ,lZ, 3t2 , H(S), A x ) such that whenever e < e , 

02, 2 [X,^] <C(O ). 
In particular, this estimate is independent of A^. 

Proof. Using the null structure equations, for each G we have an equation of the 

type 

where ip G {tr%, tr%, r), T)}, ijjjj G {x, an d V'h £ {X;^}- We also use the null structure 
equations commuted with angular derivatives: 

V 4 VV a = ^ V il ^ i2 V i3 (p + Vr Z )+ ^ V il ^ 2 V i3 (V' + ^)V i4 (^ + ^). 

U+i2+«3=« U+i2+i3+*4=« 



From the proof of Proposition 24 , we have 

II J2 ^ i2 ^ 3 p\\L^L k ms) < c(o ,n(S),A 1 ,n)ei 



il+«2+«3<2 



and 



and 



and 



V^V^vIIl-l^s) < C(O ,K{S),A l ,6 3)2 ,K)el 



«l+«2+«3<2 



^H^^M^lIl^s) < CiOo^iSl^Ajei 



U+*2+«3+*4<2 

It remains to estimate 



V^^V^VhV* 4 ^ V n V i2 V i3 ^V i4 ^. 

il+«2+«3+«4<2 ji+«2+*3+M<2 
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For the first term, as in the proof of Proposition 17, we first fix u and bound the L^ 1 L 2 (S') 
norm for each fixed u. 

II J2 v Ji # 2 v i3 ^v J4 ^|u iL2(5) 

«l+«2+*3+*4<2 

<Ce5(l + \\i)\\ L ™ L °°(s)) 2 (\\i>H\\L™L°°(s))Q2 ll v ^MU| L 2 (5) ) 

i<2 

+ Ce»(\\Vi>H\\L^L*(S)){\\^M\l%L*{S)) 

+ Ce^(\\V^\\ L ^ L 2 {S ))(\\lpH\\L^L^(S))(\\ , ipH\\LlL^(S)) 

+ Ce^ (1 + I \i/>\\LgL<*>(S))(%2 1 1 Vi3 ^l \l^lz(s)){\ \*Ph\ U|L°°(S)) 

i<2 

<C'(Oo,^(5))e^(^||V^||L r L 2 (s))(l + ||V 2 ^IU i L 2 (5)) 

i<2 



<C(0 o ,ft(SO)e5(X; l|VVslUf=L 2 (5))(l + 

i<2 

Finally, we have the term V* 1 ^* 2 V^V'ffV 14 '^. As before, we have, for each fixed 



U+i2+l3+l4<2 

<c(o ) e l (]T ||vv s |U g . La(S) )(i + \\v 2 ^\\ LiL 2 (s) ) 



i<2 



<C(O )eHj2 I|VVh|I^^(5))(1 + e^Ax). 



i<2 



Putting all these together, and using Proposition [6j we have, for each u 
||V 2 (x, u)\\ L ^ L 2 (s) 

<C(O ) + C(O , K(S), A 1; 03,2, n)e^ + C(O , K(S), A 1 )el 1 1 VVh||i^(s); 

i<2 



Taking the L 2 norm in u and using HVVtfll < C(Oo, H(S), Ai), we get 



i<2 

I|V 2 (x^)||l S l^ ( 5) <<?(Oo) +C(O ,K(S), A 1 ,0 3)2 ,^)ei 

The conclusion follows from choosing e sufficiently small depending on C , T^(S), A 1; ©3,2, 72. 

□ 

We now prove the estimates for V 2 ipu where ipH G {x? 1 ^}- 
Proposition 26. Assume 

72 < oo, 3 , 2 < oo, 72(5) < oo. 
TTien t/iere exists eo = eo(Co,72., ©3,2, 72.(5'), Ai) swc/i i/icrf whenever e < eo, 

^2,2 [*,"]< C(O ,K(S))- 
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In particular, this estimate is independent of A±. 

Proof. Consider the following equations for ipu £ {x,u}: 

v 3 ^h = p + v?7 + ipi) + ipHip + i>H(tex + i>H_)- 

As before, we commute the equations with angular derivatives: 
V 3 V 2 V^= V il ^ i2 V i3 (p + Vrj) 

«l+«2+«3=2 

+ V n ^ 2 (V i3 ^V i4 (trx + i/>h) + V^V i4 (^ + iI>h))- 

We first consider the term involving the curvature component p: 

II v^ 2 v^|Ui L2(s) 

^ + E E ii v ^ii W))(E iiv <s Piu-LiL a(S )) 

il<li 2 <2 i 3 <2 

<C(O ,^),^)J]||VV||l^(5) 

<C(O ,K(S),K). 
The term containing V 3 r? can be estimated in a similar fashion: 

|| J- V^ 2 V i3+1 77|Ui L2(S ) 
n+«2+«3<2 

^ + E E ii^ii W))(E iiv 13+ Mi^l^ ( 5)) 

U<1«2<2 i 3 <2 

<c(o ,n(s),n) W^vhims) 

i<3 

<c(o ,n(S),n,d 3 M^i)- 

We now move to lower order terms. First, we control the terms in which both ip H and 
appear: 

II Yl vV 2 v i3 ^v^ff|Ui La(s) 

*l+*2+*3+*4<2 

<C(E ll^lliV^(5))(ll^||L^(5))(E HV^IUiL 2 (5)) 
ii<2 i 2 <2 

+ ^(E I^IIW(5))(E l|V l2 ^||^L 4 (5))(E l|V^tf|UiL<(S)) 
ii<2 « 2 <1 «3<1 

+ c / (||V 2 ^ 
Jo 

+ ( ^(l|VV'||L^L2 (5) )(||^|| L ^ L ^ (5) )(||^|| i i L ^ (s) ) 

<C7(O ,^(5))(^||VVh||l^(5)) + C7 r(||V 2 ^|| i2( 5 tl ,j)(||^||^(5 I1 ,j)^ / , 

,vi Jo 
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by Propositions [23] and 25 The term with ipn and tr% can be bounded in a similar fashion: 

|| V Jl ^ 2 V 43 ^V J4 trx|UiL 2 (5) 



H+i2+«3+M<2 



<C(O ,n(S))(£\\V^H\\L~LHs)) + C l\\\v 2 MW {Su , M\^x\\l^s u , u) )d^ , 

i<l J ° 



by Proposition [23] [24] and [27| Then, we estimate the term with ip and ipH_. 

il+?2+«3+M<2 

<c(E IHI l-l- ( s))(IIV'IIlsp£-(5))(X; IIv 42 ^IIlil 2{ 5)) 

il<2 i 2 <2 



n<2 « 2 <l 
+ C(||V 2 ^||l 2 (5. u ,,j)(||^||lil-(5)) 

<c(Oo,n(s),n,Ax), 



«3<1 



by Propositions 23, 24, 25 and the bootstrap assumption (A2). The term with only ip can 
also be controlled similarly: 

|| £ V^V^V^IUi^s) 

«l+*2+«3+*4<2 

ii<2 

+^(EiHiW(5))(En vi 

i x <2 i 2 <l 

<C{0 Qi K{S),K,&i), 



il<2 

Hl»L«(S) 



i 3 <i 



by Propositions 23, 24 and the bootstrap assumption (A2). Therefore, by Proposition [6] for 
fixed m, m, 



< 



l|V 2 ^|| L 2 (5u ^ 

clival U 2 (s ,m) + c( o, ft, o^M, ao + c(o , ft(S))QT 1 1 vVhIU-^ 



(S), 



Kl 



+ C / (||V 2 ^|| L 2 (Su , J )(||(trx,^)||^(5 tl ,j)^ / . 

JO 

By Gronwall's inequality, 

IIVVHlll^(Sua) 



< ^c\\v 2 iJ H \\ms ^ + c(o Q ,n(S), n,6 3 M AO + c(o ,^(5))(g||vv^|| L ^( 5 ))j 



x exp(||(trx,V^IUi£°°(s)} 
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By Propositions 24 and 25 , the norm inside the exponential function is bounded by C(Oq, 1Z(S)). 



Thus, we have, for each fixed u 

72 



||V ^h\\l^l 2 {s) 

<C(O 0l K(S)) (\\V 2 ^ H \\ LHSo ^ + C(O , K(S),K, 3 ,2 M,Ai) 



+C(O ,Tl(S))(52\\V i iJ>H\\L° 



'L\S), 



2 



We can now take the L norm in u to get 

||V 2 ^||l2L-L2(5) 



<C(O , K(S)) [\ | V 2 M \lIl* { s 0a ) + C(O , K(S),K, O^fo], &i)e* 

+C(O MS))(J2\\^h\\lil~lhs)) 

<c(o , n{S)){\ + c(o , n{S),n, 6 3>2 , a^), 



using Proposition 23 Choosing e sufficiently small, we have 

\W 2 M\lIl^l^s) < C(O ,K(S)). 



□ 



We now estimate V 2 tr%: 
Proposition 27. Assume 

7i<00, 3 , 2 < 00, K(S) < 00. 

T/ien there exists eo = eo(C , o,'^-, 03,2,^(5)^1) such that whenever e < eo, 

O 2t2 [trx]<C(O ,K(S)). 
In particular, this estimate is independent of Ax. 
Proof. Using the null structure equations, we have 

V 4 trx = ilnf) + ip H ^ + XX- 
We also have the null structure equations commuted with angular derivatives: 

V 4 V i trx = v n # 2 v i3 v>v i4 (^ + ?M + vH i2 ^ h xS iA x- 

By Proposition 6, in order to estimate llVHrxll^i^z^s), it suffices to estimate the initial 
data and the 1 1 ■\\l%>l\ 1 l' 2 (s) of the right hand side. From the proof of Proposition 

|| Yl V^ 2 V j W 4 ^||l~l^ ( s) < C(O ,TZ(S),TZ,A 1 )e^. 

H+«2+i3+«4<2 

and 

|| Yl V^^V^VV^VlU-L^ts) < C(O ,K(S),K, AOei 

H+«2+i3+*4<2 



24 



we have 
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The only new term compared which did not appear in the proof of Proposition 24 is the 
term involving XX'- 

|| V-V i2 V i3 xV-xii^L^ ( s) 

il+«2+«3+*4<2 

< C (J2 IMIlVl£>L~(S))E HV i2 xl|L-L|L2(S))(llxl|L-L|L-(5)) 
h<2 i 2 <2 

+ c(E IHI W*-W(£ llv i2 xlU st >L i L^ ) )(l|Vxl| i5f 

ii<2 «2<1 

+ H Vil ^llL-L^L 2 (5))(llxlli2 Ls = L ^ (5) ) 
il<l 

<C(0 o ,ft(<S)), 



using Propositions 18, 19 23 and 26 The conclusion thus follows from the above bounds 
and Proposition [6] by choosing e appropriately small. □ 

We now prove the L 2 (S) control for V 2 f], thus obtaining all the 02,2 estimates. As in 
Proposition 20 where the L A (S) estimates for Vr? were derived, we need to integrate in 



the V3 direction and will not be able to gain a smallness constant. In order to get a 
bound independent of 1Z, instead of controlling V 2 r] directly, we first estimate V/x, where 
p = — div rj — p. We then obtain the desired bounds for V 2 r/ by elliptic estimates. This 
allows us to take only one derivative of the curvature components which can be controlled 
using the 7Z(S) norm instead of the TZ norm. 

Proposition 28. Assume 

K < 00, d 3 ,2 < 00, K(S) < 00. 
Then there exists e = e (O ,7l, 7Z(S), Ai) such that whenever e < e , 

2>2 [rj\<C{0 Q ,K{S)). 
In particular, this estimate is independent of A\. 
Proof. Recall that 

p = —div Tj — p. 

We have the following equation: 

V 3 p = ipifi, a, §) + ^V(V> + iPh) + ipHyi> + ip4>(trx + 4>h) + ^XX, 

The mass aspect function p is constructed so that there is no first derivative of curvature 
components in the equation. Moreover, the equation does not contain ipu- This cannot be 
derived from signature considerations alone, but requires the exact form of the V3/? equation 
as shown in Section 12.51 

After commuting with angular derivatives and substituting the Codazzi equation 



0= V^V^trx + Vir) 



«1+«2 = 1 



we get 
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^ i2 V i3 ^V i4 (trx + ^) 
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v 3 v/i= ^ i v i2 v>v i3 (p,<7)+ vii 

1 ii+«2+«3+M=2 



+ V n ^ 2 V i3 xV i4 x. 

«l+«2+«3+«4=2 

Fix m. We now estimate each of the terms. Firstly, the term with curvature: 

|| rV* 2 VV i3 (A*)lk^ ( s) 

^ + E E HV^II^ i2(S ))(E ll^ 3 (A*)iU^(S)) 
H<li 2 <2 i 3 <l 

<C(O ,K(S)) 

by Propositions 19, [23] and the definition of 1Z(S). We then estimate the nonlinear Ricci 
coefficient terms with at most one ipn'- 

1 1 Yl Vn ^ 2 V J ¥V i4 (trx + Vh) I |lil2 (s) 

U+i2+*3+*4<2 

<c( 1 +EEii v ^iit^(5 ) )(Eii vi fe tr ^ii^)) 

H<lt 2 <2 i<l 

+ C 1 1 ^ 1 1 (5) 1 1 V 2 (V>ff, trx) | \ L i L 2 {s) 

+ C\ | V 2 (r/, trx, trx) I U~£ 2 (s) 1 1 fe, trx) I UiL-(S) 

+ c / IIV^H^^, jlK^trx)!!^^, 
./o 

<C7(O , + C7 f 1 1 V\\ \ LHSul j 1 1 trx) | |l~(s u , jdu' 

>/ o 

by Propositions [18J [19J [23j 24 and 25 We control the nonlinear term with two x : 

H ^VW^xIkL^) 

n+«2+«3+*4<l 

^ + E E Hv n ^H W))(E Hv <8 xll^(5))(llxlUaL«- (s) ) 



<1 i 2 <2 



%1<L 12 



«3<1 



<C7(O ,^(5)) 

by Propositions 18, 19 and 27 Therefore, by Proposition [6j we have, for every u, 

||V//|| L oo£2(S) 

<C(O , K(S)) +C T 1 1 V 2 77| | L2(5v j 1 1 trx) IU«(S U , u )du'. 
Jo 

We now use the div-curl system 

div r) — — [A — p, 
curl rj = — a 



(35) 
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JONATHAN LUK AND IGOR RODNIANSKI 



together with the elliptic estimates in Proposition 14 to get the bound 

||V 2 ?7|| L ^ L oo L 2 (5) 



ii<l 
X 



'J2 1 1 V * 2 ^l \LgL~i?(S) + 1 1 V 13 (p, (T) \\ L ~L~L*(S) + 1 M Il^l-l 2 



12<1 



i 3 <l 



Since the 7Z(S) norm controls V*lf and V*(p, a) in L™L™L 2 (S) for z < 1, we have 
l|V 2 r7|| L£ > L oo L 2 (5) 
<C(0 o ,ft(S))(l + ^||V* 2 



i 2 <l 



This, together with (35) and Propositions 19, 23, implies that 



||V/z|| L «.La(s) 

<C(O , K(S)) + C T \\VfA\ L2{Sul j 1 1 (ipH_, trx) |Uoo (fi r jdu'. 
Jo 

By Gronwall's inequality, we have 

||V//|| LS c L 2 (S) 

<C(O ,K{S)) exp( [ U 1 1 fe, trx) IU-(S„, J^')- 
Jo 



By Propositions 23 



Thus 



By (36), this implies 



exp( / \\^H,trx)\\L°°(s u ,jdu') < C(O ,K(S)). 

\\V[4\l~l* { S)<C(O ,TI(S)). 
\V\\\ WL2[s) <C(O ,K(S)), 



as claimed. 



(36) 



□ 



Using Sobolev embedding given by Theorem [8l we improve the estimates in Proposition 



Proposition 29. 



<c{o ,n(s)). 



i<l 



Putting all the estimates in this subsection together, we obtain 
Proposition 30. Assume 

K < 00, (9 3 ,2 < 00, K(S) < 00. 
T/ien there exists e = e (C , o ; '^-5 ^3,2, 7^(S)) such that whenever e < e 0; 

o 2 ,2<c(o ,n(S)). 
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Proof. Let 

A 1 »l7(O ,^(S)), 

where C(Oq, TZ(S)) is taken to be the maximum of the bounds in Propositions [24| [25} [26} [27| 
and[28 Then Propositions 24 25 26 , 27 and 28 together show that the bootstrap assumption 



(A2) can be improved under appropriate choice of e. Since the choice of Ai depends only on 
Oq and TZ(S), we conclude that e can be chosen to depend only on Oo,1Z, C^, ^(S 1 )- □ 

5.3. L P (S) Estimates for Curvature Components. In this subsection, we estimate the 
7Z(S) norm. For this purpose, we make the bootstrap assumption 

K(S) < A 2 . (A3) 

We first prove the bounds on j3. 
Proposition 31. Assume 

1Z < oo, (9 3j2 < oo. 
Then there exists e = eo(Co, ^3,2, A 2 ) such that whenever e < e 0; 

X)ii v ^ik^ 3 (5)<c(^o). 

i<i 

Proof. Recall the V4 Bianchi equation for /3 

V,P= Y ^ 1 V i2 (p,a)+ ^V i2 (^ + ^)V l3 (^ + ^). 

U+i2=l il+«2+«3 = l 

From this we get the estimates for (3 in LlL^L 3 (S). To see this, by Proposition |6j we need 
to estimate 

|| J2 ^V l2 (p,^)+ Y ^ 1 K^ + ^h)^^ + ^h)\\lililhs). 

il+«2 = l «l+«2+*3 = l 

We have, by Proposition [THJ 

11 y ^v i2 (p^)iu^3 ( 5) 

41+12 = 1 

< e 5(i + IIVIIl-l^l-^))!^!!^^'^)!!^^ 2 ^)) 

i<2 

<C(O ,A 2 )e*K 



and by Propositions [L8j [23] and |29j 

|| Y ^ V l2 (^ + ^)V i3 (^ + V>s)lk^£3(5) 

*l+«2+i3 = l 

<C^e^(l + ||^|| i -z, r z,-(s))(^ ll Vil (^ + ll Vi2 ^ + ^)IUa^^ 4 (^)' 

il<l 12<1 

<C(O ,A 2 )el 
Therefore, 

I \§\\LiLgLS(S) < C(Ko) + C(O , A 2 , TZ)^ , 

which implies 

||^|U^3( S )<C(^o) 
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for e sufficiently small depending on Oo, A 2 and 71. We now estimate V/3. Commuting the 
V4/3 equation with angular derivatives, we have 

V 4 V/3= V' il V i2 (p,^)+ Yl V n V i2 V i3 (^ + ^)V i4 (V' + ^)- 

81+12=2 «l+i2+*3+«4=2 

After taking the L\ norm, Proposition [To] implies that 

LiLiL^s) + e« 1 1 V 4 V/3| \ L 2 L 2 uL 2 (s) ) . (37) 

The initial data is bounded by the initial data norm 

\\V§\\lZL*(S Ul0 ) <^o- 
Then, we note that by the definition of the norm 71, 

\\^ 2 §\\LiLims)<cn. 

We estimate each term in the right of side of the equation for V4V/3 in L 2 L^L 2 (S). First, 
we control the curvature term: 

|| £ V 4l ^ 2 V 43 (p,a)|| LSLiL2(5) 
n+i2<2 

ii<2 i 2 <2 
+ C{\ I I \l°°L%L±{S)) ( 1 1 (P, O-) 1 1 LlLl^(S) ) 

<c(o , A 2 )n. 



Then we bound the term with and Using Propositions 23 , 29 and 30 

|| V 1 ^V^hV 1 ^h\\lililhs) 

U+«2+«3<2 

<C(E I^IIW^(S))(E H Vi2 ^H^^(5))(E 11^*11x^(5)) 
ii<2 i 2 <2 i 3 <l 

+^(Ei^iiW^w)(Eii^^iiw^w)(Eii^Hiiie.iiL"(fl)) 

il<2 i 2 <l i3<2 

+ c(||v^lU^ (g) )(||v(^,VH)IU|^(s)) 

<C(O ,A 2 ,ft). 

Since ^ satisfies stronger estimates than either or ipH_, we have 

II ]T v^ 2 v 43 (^ + v*)v i4 (v> + ^)lk^(s) 

«l+*2+«3<2 

<c(o , a 2 ,tz). 

Putting the bounds together, we have 

\\^4V§\\ LiLlL 2 iS) <c{o ,A 2 ,n). 

Thus, (37) implies that 

1 1 V£| \lil~lhs) <C (n + e*K + e*C(O , A 2 , K)) . 
The proposition follows from choosing e sufficiently small depending on G , 7Z, 0^ 2 , A 2 . □ 
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Once we prove the estimate of lZ(S)[fi] independent of the 1Z norm, we get the following 
improved bounds on the Ricci coefficients: 

Proposition 32. Assume 

1Z < oo, 3>2 < oo. 
Then there exists e = eo(Co, TZq, &3,2, A2) such that whenever e < €q, 

Oo,ooW<C(0o,^o), 
ll2 y <C{0 ,Ko), 

Y,OiAx,u,trx] <c{o ,n ). 



i<l 



Proof. This follows from substituting the bound for TZ(S)[(3] from Proposition 



estimates from Propositions 18, 19, 21] and 22 
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into the 

□ 



Using this improvement, we prove the 7Z(S) estimates for p and a. 
Proposition 33. Assume 

1Z < oo, 3j2 < oo. 
Then there exists e = e (O , 1Z , 71, 3:2 , A 2 ) such that whenever e < e , 

E II V*(p, a) \\l°?l%>l 2 (S) < C(O , n ). 

i<l 

Proof. Consider the V4 equations for p and a: 

V 4 (p, a) =V/3 + v>(A P, o) + V^v(^ + M + ^ii) + M + i>xx- 

After commuting with angular derivatives, we get 

V 4 V(p,cr)= V il # 2 V i3 (/3,p,o-)+ E V^V^V^ + W 

+ E i; ll V i2 ^V l3 xV i4 x- 

U+i2+i3=l 

By Proposition p| in order to estimate V J (p, a) in L^ L^ L 2 (S'), we need to estimate 
V iV t {p ) a) in L^L^L 2 (S). The term with curvature can be estimated by 

|| V 4l # 2 V i3 (/3,p,a)|| Lii2(5) 

ii+«2<2 

<C^ (E 1 M Il-l- ( 5)) (E 1 1 Vl2 ^ A *) 1 1 *W>) 

ii<2 i 2 <2 

+ Ce5(||V^|U^ (5) )(||(/3,p,a)|U ii 4 (s) ) 
<C(O ,K(S),1l)e^K 
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by Proposition 23 The nonlinear Ricci coefficient term with at most one ipn can be controlled 
by 

II V il ^ 2 V i ^V i4 (^ + ^)|U i L 2 (5) 

<e* 1 1 Yl V 11 ^ 2 V^V i4 (if, + M\\li^ { s) 



ii<3 i 2 <2 
i 



+ Ce2(||V^|U^4 (5) )(||V(^,V H )l|L|^(5)) 
<C(O ,^(5),^)e5 



by Propositions 23 and pOj The remaining term containing \\ can be estimated using 



Proposition 32 



^V^V^V^xIIl^os) 

U+«2+*3<l 



< 



\x\ 



(S), 



il<2 



i 2 <l 



+ C| |^l U^L-(S) 1 1 V^| |l-L2(S)(| |^| |l|L-(5)) 2 

<c(o ,n ). 

Therefore, by Proposition [6] 

1 1 V*(p, *) I U 2 ^,*) <C(^o, ^ ) + ^C(O , Ko, K(S),K). 



i<l 



By the bootstrap assumption (A3) on 1Z(S), we can choose e small depending on O , 1Z , 
TZ, and A 2 to conclude the Proposition. □ 

Finally, we prove the bounds for the Gauss curvature. This will be used in the next 
subsection to carry out elliptic estimates. 

Proposition 34. Assume 

1Z < OO, 03 ; 2 < OO. 

Then there exists e = eo(Co, ^0, @3,2, A2) such that whenever e < €q, 



i<l 



Proof. Consider the Gauss equation: 



K = —p + ipip. 



We estimate each term on the right hand side. By Proposition [33 

\l™l™l 2 (s) < C(O ,TZ ). 



By Propositions 23 and 32 



II^IIl-l^(s) < 11^11^-1,4(5) < c(o ,n ). 



□ 
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We can thus close the bootstrap assumption (A3) to prove the following estimates for 
7Z(S), under the assumption that 7Z and (9 3)2 are bounded. 

Proposition 35. Assume 

1Z < oo, 3>2 < oo. 
Then there exists eo = eo(Co, 7Zo, TZ, ^3,2) such that whenever e < e , 

K(S)<C(O ,K ). 

Proof. Let 

A 2 >L7(0 ,^o), 

where C(Oq, TZq) is taken to be the maximum of the bounds in Propositions 
Hence, the choice of A 2 depends only on Oq and TZo. Thus, by Propositions 
the bootstrap assumption (A3) can be improved by choosing e sufficiently small depending 
on O ,TZ ,TZ and <f) 3)2 . □ 

Using Proposition [35l we improve our estimates on the Ricci coefficients in Propositions 
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33 


and 


34 


31 




33 


and 


34 



23 and 30 to get the following: 



Proposition 36. Assume 

1Z < 00, 3j2 < 00. 
Then there exists e = e (O ,7l ,7l, 3 , 2 ) such that whenever e < e , 



i<2 



5.4. Elliptic Estimates for Third Derivatives of the Ricci Coefficients. We now 

estimate the third angular derivatives of the Ricci coefficients. Introduce the bootstrap 
assumption: 

63,2 < A 3 . (A4) 
The bounds for the third derivative of the Ricci coefficients cannot be achieved by the 
transport equations alone since there will be a loss of derivatives. We can however combine 
the transport equation bounds with the estimates derived from the Hodge systems as in [18] , 
[8], [IH]. We first derive the control for some chosen combination of V 3 (ip, ipH, 4>h) + (/?, p, a) 
by the transport equations. Then we show that the estimates for the third derivatives of 
all the Ricci coefficients can be proved via elliptic estimates. We begin with the bounds for 
V 3 tr% and V 3 %: 

Proposition 37. Assume 

11 < 00. 

Then there exists e = eo(Co> 1Z , 1Z, A 3 ) such that whenever e < e , 

6 3t2 [trx,x}<C(O )(l + 1Z[P}). 

Proof. Consider the following equation: 

V 4 trx = XX + trxO + M, 
After commuting with angular derivatives, we have 

V 4 V 3 trx= V il # 2 V i3 trxV i4 (^ + ^) + V* 1 ^ 2 V 43 xV i4 x- 

U+i2+«3+«4=3 il +i2+23+«4=3 
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We estimate term by term. First, we bound the term with tr% and ipH- Integrating in the 
u direction, applying Sobolev Embedding in Propositions [8] and [9] and using Proposition 36 

we get 

|| £ VV 2 V* 3 tr X V^lk^ (s) 

n+i2+l3+M<3 



< C ^C£ E ii v ^i!W))£ nv i3 tr X ii^^(5))(E iiv i4 ^iuil 2( 5)) 

+ liV ll V'iliV^(5))(E HV l3 tr X || LrL2(5) )(^ ||V^MkzW 



ii<2 « 2 <3 

<C(O ,^ ) e i(l + A 3 ). 



i3<3 



i 4 <l 



Since ?/> satisfies stronger estimates than ipn, we have the same bounds for the term with 
tr% and if): 

n+«2+*3+*4<3 

<C(Oo,^ )e5(l + A 3 ). 
Finally, we consider the term with XX'- 

|| £ VV 2 V i3 xV J4 x|U lL2( 5) 

«l+i2+«3+*4<3 



<^(£En vJ1 

ii<2 i 2 <3 



:Eii v ^ii^(^))(Eii v ^ii^(^))^ 

23<2 24<3 

<C(Oo,^ )(l+ /^||V ^ x||L 2 (5 ^I , a ,))(l|V 3 Xl|L 2 (5^))^ , )• 



(38) 



i<2 



We now use the Codazzi equation 



and apply elliptic estimates in Proposition [15] to get 

l|v 3 xlU 2 (5) 

<C(0o,ft o )(E \\™X\W {S) + 2 11^/311^(5) 



i<3 



i<2 



(39) 



+ ^ \\^ipV l2 (iP + ip H )\\ms) + \\x\\LHs)) 

ii+«2<2 



Notice that we can apply elliptic estimates using Proposition [15] since we have bounds for 
the Gauss curvature by Proposition 34 Therefore, 

\Eii v ^ii^(^))(iiv 3 xiu 2( 5„, a ,))^ 

i<2 

<c(o ,n )(i+ AX)llv < ^||^^,)(||v%rxll^.o)^ , + ^)- 

i<2 
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Gathering all the estimates, we get 
l|V 3 trx|U 2( s U]a ) 

<c(o ,n )(i + e^A 3 + r(X;ilv i2 xlU 2 (5„ u ,))(ll v3tr xllL 2 (5„ tt ,))^ + W])- 

J ° i 2 <2 

Gronwall's inequality implies that 
l|V 3 trx||L2 ( s^) 

<C(O ,n )(l + e^A 3 + K[P])eM /~(E liV'xIU^,))^') 
<C(O ,n )(l + e^A 3 + ft[/?])exp(es ^11^x11^^(5))) 

i<2 

<C(O ,7l )(l + e^A 3 + n[/3]). 
By choosing e sufficiently small depending on Oo,TZ and A 3 , 

l|V 3 tr X || L2( ^ ) <L7(O ,^ )(l+^[/3]). 
This, together with (39), implies that 

\\v 3 x\\ L ^Lims) <c(o )(i + n[ 



□ 

We now prove estimates for V 3 //. To do so, we first prove bounds for second derivatives 
of p = — div 1] — p and recover the control for V 3 rj via elliptic estimates. 

Proposition 38. Assume 

n<oo. 

Then there exists e = e (O , 1Z , 71, A 3 ) such that whenever e < e 0; 

3j2 [p,r]) <C{O )(l + e^d 3>2 + K). 

Proof. Recall that 

p = —div Tj — p. 

and satisfies the following equations: 

V 4 /i = V(A P, cr) + VV(^ + *p H ) + + + M + i>XX, 

It is important to note that (3, ipn_ are absent in this equation. This cannot be derived from 
signature considerations alone, but requires an exact cancellation in the equation for V 3 p as 
indicated in Section 12.51 

After commuting with angular derivatives, and substituting the Codazzi equation 

/3 = V X + ^ n V i2 ^, 

we get 

v 4 v 2 /i= E v n ^ 2 v i3 ^v 44 (p,a) + E v^v^v^ + ^M 

*l+«2+*3+*4=2 il + t 2 +j3+*4=3 

il+l2+«3+i4=3 
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The term with curvature can be estimated using Proposition 36 by 

II ]T v ii v i2 v i3 VV^(p,^)|U» z i i2(S) 

U+«2+«3+«4<2 
?i<2i 2 <3 i 3 <2 

<C(0 O , n )e^ l|V*(p, a) 1 1 



j<2 

1 



<C(Oo,^ )e5^ 



We next consider the term with two x's. By (38) in the proof of Proposition 37, we have 
|| VV 2 V i3 xV i4 xl \l~i> & i?(S) <C(Oo,^q)(1 + o 3 , 2 [x]) 

il+«2+«3+M<3 

Applying Proposition [37] we get 

II ^ v li ^ 2 v l3 xV i4 xllL-L^(5) < c(o ,^o)(i + n). 

ii+?2+*3+«4<3 

We then move to the remaining terms with at most one ipH- First, we look at the terms that 
do not contain ipuV^ip. These are the terms 

«l+«2+«3+M=3,l3<2 

and 

«l+«2+«3+*4=3 

The first term can be estimated using Proposition [36] by 

|| ^H l2 ^ l3 ^ U M\L^L k LHS) 

«l+*2+«3+«4=3,i3<2 

^ e *(E E ii v ^n WWE iiv^Hl-l-l4(5))(E \\^^h)\\l^( S) ) 

ii<li 2 <3 «3<1 i4<3 

+^(Eii^iiW i -w^Eii^ a ^ii i ^ ia w)ii^'^)iiw i -(^ 

ii<3 i 2 <2 

<C(0 ,^o)e5(l + A 3 ). 



The second term can be controlled using Proposition [36] by 

|| £ VV 2 V^V i4 V>ii^^(s) 

*l+«2+«3+i4=3,i3<2 

< Gel E E n^n WWE hv^iitol 4( 5))(E iiv^iil- 

jl<l«2<3 *3<1 «4<3 

<C(Oo,^ )e5(l + A 3 ). 

We now bound the terms ipH^ip- If ^ € {^ T Xi ^ r x}) we can estimate in a similar fashion as 
before, since we have L^L^L 2 ^) estimates for V 3 (tr\;, trx): 

| |^V 3 (tr X , trx) Wl^lW(s) <C(O , TZ )e^A 3 . 



NONLINEAR INTERACTION OF IMPULSIVE GRAVITATIONAL WAVES 65 

The remaining terms are of the form ipH^ 3 {VyV)- The difficulty in estimating these terms 
is the fact that using the (D 3 ^ norm, V 3 r? and V 3 ?7 can only be estimated in L 2 (H) but not 
L 2 (S). Thus we need to estimate both V 3 (rj,r]) and ipn in L\ and will not have an extra 
smallness constant . Therefore, instead of bounding V 3 (r], r\) with the norm, we apply 
elliptic estimates and control V 2 /i in L^ L^°L 2 (5'). 
By the div-curl systems 

div r] = —fi — p, curl 7] = a, 
div rj = p + p, curl r] = —a, 



and the elliptic estimates given by Propositions [14] and 34 , we have 

\\V 3 V \\ms) < C(J2 IIVVII^(S) + ^ ||V J (p,a)|| L2(s) + WvWms)), (40) 

i<2 i<2 

\\v 3 v\\ms) < c(J2 llvVlU2 (S ) + W^&^Wms) + ll!zlU 2 (5))- 

i<2 i<2 

This implies 

\\^HV 3 {ri,V)\\LlL2(S) 
<C\ \ll> H I \l%L<*>{S) 1 1 V 3 (r/, 77) I |l|L2(5) 

i<2 i<2 

<c(o ,n Q )(i + e^A 3 + n). 

Hence, gathering all the above estimates, by Proposition [6j we have 

||V 3 /i|| LS o LrL2(5) < c(o ,n )(i + e^A 3 + n). 

By choosing e sufficiently small depending on A3, we have 

||v 3 /i|| LS o LS o L2(s) <c(o ,n )(i + n). 



Therefore, by (40), we have 

\\V 3 V \\ L? o L 2^2 is) < C(O ,K ){l+K), 

and 

\\V 3 v\\L~LiL*( S )<C(O ,K )(l+K). 

We now estimate V 3 a;: 
Proposition 39. Assume 

n<oo. 

Then there exists e = e (O , 1Z , 71, A 3 ) such that whenever e < e 0; 

3 , 2 [k,uj_,J] < C(O ,7l )(l + 7l[P)). 



□ 
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Proof. Recall that u} is defined to be the solution to 

„ + 1 



with zero initial data and 



K := -Vw+* Vw 1 " - -/?. 



By the definition of u;t, it is easy to see that using Proposition [6 

X)ll V V|U^^(s)<C(0o,^o) 



j<2 



for sufficiently small e. In other words, u} satisfies the same bounds as ipn- in the remainder 
of the proof of this Proposition, we therefore also use ipH to denote ufl . 
Consider the following equation for k: 

V a k = i)(p + a) + (i) + i) H )V{i} + ^h) + + ^)V(V> + M + ^(ip + ^h){^ + ^h). 
After commuting with angular derivatives, we get 

V 4 V 2 /i= y h ^ i2 V i3 ^V i4 (p + a)+ V il # 2 V i3 (^ + ^)V i4 (^ + ^)- 

11 +12+13+14=2 U+l2+i3+M=3 

We estimate k in L^L^ L? (S) . Fix u. By Proposition |6j for each u we need to estimate 
the right hand side in L\L 2 (S). After taking the L 2 norm in u, we thus need to control the 
right hand side in L^L^L 2 ^). The term involving curvature has already been estimated in 
Proposition 38 and can be controlled by 

«l+«2+i3+«4=2 

Thus, 

|| Yl V il ^ 2 V i3 ^V l4 (p + <j)\\LiL iL z(S) < C(0 ,Ko)e*K. 

For the other terms, it suffices to consider 

E V H ^ l2 V i3 ^ H V u ^H 

*l+«2+«3+*4=3 

since ip satisfies stronger estimates that either ipjj or ipH- To this end, we have 
|| Y V^V^V^k^s) 

11+12+13 +«4<3 

<Cei 1 1 ^ 12 ^h^h\\lIlIlHs) 

n+«2+«3+*4<3 

< C ^(£ E ll V ^li W*( S )XE li^ 3 ^iU»L|L 2 (5))(E W^mWlIL^LHS)) 
ii<2i2<3 *3<3 «4<1 

+ ^E E ll V ^HW 2 W )( E H^ 3 ^|U|^^(5))(E HV^II^L S L 2 (5)) 
ii<2i2<3 *3<1 «4<3 

<C(Oo,^o)e^(l + A 3 ) 
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by Propositions 36 Therefore, by Proposition |6j 

\\VW\lilzli{s) <C(O ,K )(l + e 1 i('R + A3)). 
By choosing e sufficiently small depending on 1Z and A 3 , we have 

II V «|| iALF&iS) < C(O ,TZ ). 
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Consider the div-curl system 



div Vw = — div k — @i 
curl Vw = 0, 



div Vw = —curl k — -curl (3, 
curl Vw = 0. 



By elliptic estimates given by Proposition [14] and [34| we have 



□ 



In the remainder of this subsection, we consider the third derivatives of the Ricci coeffi- 
cients that are estimated by integrating in the u direction. We need to use the fact that the 

are independent of A3. We now estimate V 3 tr% 



37 



38 



39 



estimates derived in Propositions 
and V 3 x: 

Proposition 40. Assume 

n<oo. 

Then there exists e = e (O ,7l ,7l, A 3 ) such that whenever e < e 0; 

3 , 2 [trx,x] < C(O ,n )(l +K\g). 
Proof. Consider the following equation: 

V 3 trx = XX + trx(trx + Vg) 
After commuting with angular derivatives, we have 

V 3 V 3 trx= V^V^V^trxV^^trx + VfL) + J2 V^V^xV* 4 *, 

«l+*2+i3+*4=3 il+l2+i3+i4=3 



Fix u. We estimate term by term. First, by Proposition 36 

|| V 4l ^ i2 V i3 trxV J4 trx| \ W s) 



«l+«2+23+M<3 

<^(EEn vS 

n<2 i 2 <3 



(J2\\V H trx\\ LH s u ,Jdu' 

1 I i 4 <3 



<C(O ,1l )(l + f \\V 3 tTx\\ms u ,jdu'). 
Jo 
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Then we bound the terms with one We separate the cases where ipH_ = an d — X- 
First, for i/>ff = ur. 

|| J2 V 4l V i2 V i3 trxV 4 V||LiL 2 (5) 

<c(E E i iv^ii^ l2(5) )(E 1 1 v i3 trxi u W lh5))(E 1 1 v >i iww) 

ii<2i2<3 «3<1 i4<3 



L2(S)> 



(E H Vi3 *^',j)(E liV-trx|U 2( 5 tt , j)^' 



<C7(0 ,72o)(l +^] + / C£ llV^adk*^, j)(Ell Vi3tr xiU 2 (V,j)^')> 



«2<1 



*3<3 



where we have used Propositions 36 and 39 Then, we consider the term with one tpn, where 
= X- 



|| J2 V 1i ^ 1 ^xV 1 %\lilhs) 

<^(E E h v ^h WW f (£ n vi3 iii^j)(E 11^11^(5^))^ 

ii<2i 2 <3 ^° i 3 <l i 4 <3 

+ <?(£ E il V ^H WW AE HV i3 trx|U,(5 u , j)(E llV**xllL>(s^,)d«' (41) 

ii<2i 2 <3 £ 3 <1 i 4 <3 

<c(o ,^o)(i+ r(Eii v iii^j)(ii v3tr ^i^(^j)^) 
+c(o ,^o) r(En vitr ^i^(^,j)(ii v3 iiu 2 (^,j)^ / > 



using Proposition 36 In order to control this, we need to use the Codazzi equation 



div x = 2 Vtr ^ + §_ + W + M 



and apply elliptic estimates using Propositions 15 and 34 to get 



l|v 3 xlU 2 (5) 

i<3 i<2 

L 2 (S) + ||x||l 2 (S)) 

ii+«2<2 



(42) 
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Using (42), we can bound the second term in (41): 



£||V i trx|U,(^,j)(||V 3 xlU 2 ( V> j)^ / 

i<l 

<C f U (J2 \W l trx\\LHs u , j)(\\Vhrx\\L> { s u ,Jdu' 

JO 



i<l 

+ C7(^ \\VHtx\\l~i H s))(£ \\^ 2 g\Lims)) 

il<l i- 2 <2 



h<l i2+«3<2 
« u 



<C(O ,K )(l+K[§} + / HVHrxll^^, j)(||V 3 trx|U 2( 5 u , j)^) 



ii<l 

This, together with ( [41] ) , implies that 

|| £ V^ 2 V* 3 trxV l4 xlkL 2 (5) 

<L7(0 ,^o)(l + ^[^]+ ^(^||V ^ (trx,x)l|L 4 (5 tl ,J)(||V 3 trx|U 2 (5 Il ,J)^ , )• 

- 70 i<i 

Finally, we estimate the term with two ipt^s. We note that the only such term is of the form 
XX- We control this term using (42): 

|| J2 ^^V^xV u x\\lilhs) 

ii+i2+h+U<3 



<^(E E ii^ii WW f(E n vi3 iii^(^,j)(E ii^m^j)^' 

ii<2« 2 <3 ^° i 3 <l i 4 <3 

<C(Oo, + C AE II 1^ j) (I I V 3 trx| |l 2(Su , 
■ / ° i<l 

+ C(E IIv 1i xIIl^( 5 ))(E liv^lk^)) 

tl<l i 2 <2 

+ c(Ell vn iH^ 4 (^))( E IIv^v 13 (^ + ^)IIl^(5) + IIaIIl^(5)) 

pu 

<c(o ,n )(i + n[§] + / (Ell v llU 4 (^j)(ll v3tr ^IU 2 (^j)^ / )- 

Therefore, by Proposition [6j we have 
l|V 3 trx||z~L2(5) 

<C(O ,TZ )(l+n[§] + /\Ell Vl ^'^)ll^(W)(ll V3tr XlU 2 (^,j)^')- 



Kl 
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By Gronwall's inequality, we have 
l|V 3 trx||L-L2(5) 
<C(0 o ,fto)(l + ft[£J)exp( r&\\V i WM\\LHs v ,Jdu') 



" u »<i 

<C(O , K )(l + K\g) exp(C(O ,n )) 
<C(0 ,Ko)(l + K\p)). 



By (42), this implies 

||V 3 x||l-l^(5) <C(O ,K )(l + K 



□ 



We now control V 3 i]. 
Proposition 41. Assume 

n<oo. 

Then there exists e = e (O , 1Z , 1Z, A 3 ) swc/i t/iat whenever e < e , 

<C7(O )(l + ^). 

Proof. Recall that 

fj, = — div rj — p. 

We have the following equations: 

v 3 // = ^(p, o-) /3) + ^v(-0 + + + t/}ip(ip + iPh) + i>xx, 

After commuting with angular derivatives, we get 

V 3 V 2 /i= V il V i2 V^V i4 (p,CT)+ v^^v^vv**^ + V's) 

n+«2+«3+*4=3 

We estimate every term in the above expression. First, we bound the term with curvature: 
|| £ V 4l # 2 V l ^V i4 (p,a)||^ iiL2(s) 

21+i2+i3+*4=2 



<c(E ii^nw^(5))(E iiv i2 (p,*)iu^ ( s)) 

ii<2 i 2 <2 

+ c(^||^||iv^ i o 0(s) )(^||v i ^|U^^ (s) )(^||v^( / 5, ^)IU r ^ (s) ) 

H<2 8 2 <1 « 3 <1 

<C(0 ,^o) ^ \\V l (p,a)\\ L ~ LlLHs) 

i<2 

<c(o ,n )n. 



We now move on to the terms with the Ricci coefficients. Notice that by Propositions 37, 38 
40l all the terms of the form V 3 ?/> except V 3 r/ have been estimated. Thus, by Propositions 
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«l+«2+i3+«4=3 

<^(EEii v ^iiW 2 (5))(Eii vi3 ^'^iiw^^)) 



il<2i 2 <3 i 3 <2 

x 1 1 V i4 (trx, try_, 1 1 (S)) 

i 4 <3 



<l7(0 ,^o)(1+^) + l7 / (EHV i (V',^)IU 4 (5 u ,jliV 3 !Z |U 2 ( Su , i j^ / . 



i<l 



We control the last term using the div-curl system 

div T] = fi + p, curl r\ = —a. 
Applying elliptic estimates using Propositions [14] and |34j we get 

||V 3 r/|| L2(s) < C(O 0l TZ )(J2 I|V>||l3(s) + E H V ^' & )W^(S) + \\v\\l*(S)) 



i<2 



i<2 



Thus, we have 



i<l 



<c(o ,n )(i + n+ / (Ellv*(^«IU«(^)l|v 2 A*|| L » (fl L /ia )d« / )- 



i<l 



Hence, by Proposition |6j we have 
l|VV||^ L2( 5)<C(Oo,^o)(l + ^ 

By Gronwall's inequality, we have 



(Eii vl ^'^)ii^(^,j)ii v vn^(^,j^ / ). 



i<l 



[4\l™l 2 (s) < C(O , 1l )(l + K) exp( / V 1 1 V<(V, ^ I U* ( s u , j I U* ( s„, jdu')- 



By Proposition 



36 



l 4 (s u ,u) is controlled by C(O ,TZ ). Therefore, 
llV^IU^s) <C{O ,n ){l + K). 



i<l 



The desired estimates for V 3 r/ thus follow from (43). 

We finally prove estimates for V 3 w. 
Proposition 42. Assume 

K<oo. 

Then there exists eo = eo(Co, T^-o, A3) swc/i i/icrf whenever e < eo, 

03, 3 [/t,w,a;t]<C?(Oo)(l + ^). 
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Proof. Recall that a/ is defined to be the solution to 

1 



with zero initial data and 



VgC^ = -<T, 



K := Vw +* Vw f - -p. 



By the definition of uj\ it is easy to see that using Proposition |6j 

J2\\^\\LiL^( S )<c{o ,no) 

i<2 

for sufficiently small e. In other words, satisfies the same estimates as ipn- in the 
remainder of the proof of this Proposition, we therefore also use ipa to denote aA Consider 
the following equations: 

V 3 k = V(P + o) + (V> + V^)V(V> + ?M + (V> + ^)V(V> + ?M + W + + iPh), 

Commuting with angular derivatives, we get 

V 3 V 2 k= J2 V n ^ 2 V i3 ^V i4 (p + (j) + E V il # 2 V i3 (^ + ^)V i4 (^ + ^). 

*l+j2+*3+*4=2 il + t2+i3+*4=3 



Fix u. The term involving curvature has already been bounded in Proposition 41 and can 
be controlled by 

|| £ V Jl ^ 2 V 43 VV J4 (p + a)|| L i L2(s) < C(0 ,^o)^. 
For the terms with only Ricci coefficients, notice that the third derivatives of all the Ricci 



coefficients except u and have been estimated. Thus using Propositions 36, 37, 38, 39 



40 and 41 , we have 



xE 1 1 v * 3 faM i iws)) (E 1 1 v * 4 ^) i 



i i e yn ^ 2 v * 3 ^ ^) yi4 ^> ^ lUi^(^) 

U+«2+i3+«4<3 

<c(EEii v ^n^x 

ii<2« 2 <3 « 3 <2 i 4 <3 

+ c(| I v 3 xl \lilhs)) (E 1 1 v ^ ^) Ik^)) 

i<2 



' i<l 

<C(0 ,^o)(l + E HV^Ik^W + I|V 3 x||l^(5) 



ii<2 



+ f ||V 3 (w, w t) | U 2(5u , j(E 1 1 V i2 (V, 1>e) \\LHs u ,jdu' 
Jo i2 <i 
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Therefore, by Proposition [6j 

I|V 2 k||l S ol2 (5) 



<C(O ,n )(l + K + WVHhWizl^s) + I|V 3 x||l^(5) 

l<2 

1 1 vV, w t) i | i2(5v a) 1 1 v* 2 ^, Va) I |^(s u ,, 



«2<1 



By the following div-curl system: 



div Vw = div k + 2 div Pi 
curl Vo; = 0, 
div VuA = curl k + -curl 0, 
curl VuA = 0, 



we have, using Propositions 14 and 34 



||V 3 (u;V)|| L2(5) < C7(||V 2 «:|| L2 (5) + 1 1 V 2 /3| | L 2 (S) + \\V(u,J)\\ L , {s) ) 
Applying this to the estimates for V 3 k and using Proposition [36] we get 



(44) 



U<2 



+ f 1 1 V 2 k| U^, j (£ 1 1 V " M « \\L H s u ,jdu>). 

J ° ~ fe<l 



By Gronwall's inequality, and using Proposition 36 

I|V 2 k||l S cl 2 (5) 



<C7(O , ^o)(l + ^ + X: 1 1 Vtyffl k^(5) + I|V 3 x||l2l2(5)) exp( / | |^ s |U~(s tt , J**') 

i<2 ^ 

<l7(Oo,^ )(1 + ft + J] ||VVh|U S l 2 (s) + ||V 3 xlU Si2(s) ). 



i<2 



Now, taking the L 2 norm in u, we get 



|V 2 /d 



<C(C )(1 + ^ft + l|VVff|U|^(S) + ||V 3 x||l|l S l 2 (5)) 



i<2 



<C(O ,^ )(l + e2^ + 7e[/9]), 



where in the last line we have used Propositions [36] and [37| By choosing e sufficiently small, 
we have 

||V 2 k|| L 2^ l2(s) <C(0 ,^o)(l + 7e[/3]). 



Therefore, by (44), we have 

||V 3 (o;,a;t)|| L ^2 i2(s) <C(Oo)(l + 7e 



□ 
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Putting these all together gives 
Proposition 43. Assume 

n<oo. 

Then there exists eo = eo(Co, T^o, TV) such that whenever e < eo, 

6 3 , 2 [trx, x, k, 00, J] < C(O , K )(l + K[/3]) 7 

03 !2 [trx,x,&^,^} <C(Oo,Ko)(l + K[ 



and 



Proof. Let 



6 3)2 [v,v,v,v] <C{O ,TZ )(l + TZ). 



A 3 >C(C ,^o)(l + 7e), 



where C(Oq, 1Zo)(1+1Z) is taken to be the maximum of the bounds in Propositions 37, 38 



40 



38 



41 



39 



[42| Hence, the choice of A 3 depends only on G , TZq and 1Z. Thus, by Propositions 



39 



37 



40, 41 42, the bootstrap assumption (A4) can be improved by choosing e sufficiently 



small depending on Oq,1Zq and 1Z. □ 



6. Estimates for Curvature 

In this section, we derive and prove the energy estimates for the curvature components 
and their first two derivatives and conclude that TZ is controlled by a constant depending 
only on the size of the initial data. By the propositions in the previous sections, this shows 
that all the O norms can be bounded by a constant depending only on the size of the initial 
data, thus proving Theorem |4| In order to derive the energy estimates, we need the following 
integration by parts formula, which can be proved by a direct computation: 

Proposition 44. Suppose <f>% and 02 ore r tensorfields, then 



0lV 4 02 + / 02V 4 01 = / 0102 - / 0102 + / (2u> - trx)4>l4>2 

D u ,u JD UA JH u (0,u) JILo(.0,u) JD u ,u 



0lV 3 02 + / 02V 3 01 = / 0102 - / 0102 + / (2w - ^)010 2 . 

D„,u JD u ,u JH u (0,u) JH (0,u) JD u ,u 

Proposition 45. Suppose we have an r tensorfield ^0 and an r — 1 tensorfield ^0. 

= - f O? + 2 ) (1) (2 V- 
Using these we derive energy estimates for p, a in L 2 (H U ) and for f3 in L 2 (.ff u ). 



NONLINEAR INTERACTION OF IMPULSIVE GRAVITATIONAL WAVES 

Proposition 46. The following L? estimates for the curvature components hold: 



75 



i<2 



<E(ii^*)iiW-j + 11^11 W.o)) 

£ vi (^))( E rv i2 (v+^)v i3 (p^))ii wl( 5) 

i<2 u+i2+J3<2 



i<2 
+ 



+ 
+ 



+ 
+ 



i<2 ii+i2+«3+*4<3 

(£ V (/>,*)) ( E ^ V 13 x V- X) 1 1 L\L\L X (S) 

i<2 ii+i2+i3+«4<2 



i<2 



H+i2+»3<2 



Ev^)( e r^ i2 ^Hp,mLiL k L H s) 

i<2 u+«2+j3<l 

E V^)( E V^V^rx + ^)V i4 (^ + Vh))|| 



i<2 ii+i2+*3+«4<3 

Proof. Consider the following schematic Bianchi equations: 

V 3 <x + div *§_ =i)b + E ^ V i2 ^V i3 ^ + ^XX 

U+i2+i3<l 

V 3 /5 + div/3=# + E ^V^V^ + V'XX' 
n+«2+«3<i 

V 4 /3 + Vp -* V<7 =ij)(p, a) + ^ V i2 (^ H + trx)V i3 (^ + trx), 

Commuting these equations with angular derivatives for i < 2, we get the equation for 
V 3 V l a, 

V 3 V i <7 + div * V l /3 

= E V' il V i2 (V' + ^)V i3 (p,<T)+ E V^^V^^V^^ + ^h) 

H+«2+»3<2 ll+l2+j3+*4<3 

+ E ^ n V i2 ^V i3 xV i4 x- 

*l+«2+«3+M<2 

Notice that in the above equation, there are terms arising from the commutator [V*,div ]. 
These can be expressed in terms of the Gauss curvature, which can be substituted by — p — 
|trxtrx and rewritten as the terms in the above expression. The equation for V 3 V*p has a 
similar structure: 

V 3 V l p + div V*£ 

= E V' n V l2 (^ + V'//)V l3 p+ E V n ^ l2 V i3 ^V i4 (^ + ^) 

«l+«2+«3<2 «l+«2+i3+*4<3 

+ E ^V^V^xV^x- 

U+i2+«3+*4<2 
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We have the following equation for V4V*/3: 
V4V73 + VV\5 -* VV*<7 
= V> il V i2 V>V* 3 Go,<7)+ ^V i2 pV i3 (p,CT) 

+ V n ^ 2 V i3 (^ + trx)V i4 (^H + trx). 

«l+«2+*3+«4<3 

As a shorthand, we denote by F\ the terms of the form 

Fx := ^ V il V i2 (V' + V'H)V i3 (p ) o-)+ v^v^v^ + V'h) 

ii+i2+i3<2 ii+«2+i3+*4<3 

+ ^ n V i2 ^V i3 xV i4 x, 

n+«2+«3+*4<2 

and by F2 the terms of the form 

F 2 := ^ 1 V i2 ^V i3 (p,a)+ ^ V^V^pV* 8 ^,*) 

»i+ia+»8<2 ii+i2+i3<l 

n+*2+*3+i4<3 

Applying Proposition 45 yields the following identity on the derivatives of the curvature. 

< p, -VV> +* VVV > 7 + < V*£, F 2 > 7 

< div V* P, V*p > 7 + < div * VV3, V'a > 7 + < V*/3, F 2 > 7 
y - < VgV'p, V'p > 7 - < VgV^, W > 7 + < V*£, F 2 > 7 + < V*(p, a), F x > 7 . 



Using Proposition 44 we have 



I <v^,v 4 v^> 7 =W |v^| 2 - / IvW + ik^-^iv^HUi^cs). 

Substituting the Codazzi equation 

£ = W>h + V# + <M 

for one of the /3's, we note that the last term 

-hr X )\V% 2 \\ LlLlL , {s) 
is of the form of one of the terms stated in the Proposition. We call such terms acceptable. 



Also by using Proposition 44, we have 



J <V\p,a),V 3 V l (p,a) > 7 
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The last term 



(w- UrxW(p,a)\ 2 \\ L{iLiL i {s) 



is also acceptable. We thus have 

/ |V*£|? + / |VW)I? 

< [ |V\0|* + / \V\p,a)\ 2 + \<V^,F 2>7 \ + \<V l (p,a),F 1 >, 



+ acceptable terms. 

We conclude the proposition by noting that the structure for F\ and F 2 implies that 



and 



are also acceptable. 



<V%F 2 > 1 \ + \<V i (p,d-),F 1 >, 
<V%F 2 > 1 \ + \<V i (p,&),F 1 >, 



□ 



To close the energy estimates, we also need to control /3 in L 2 (H) and (p, a) in L 2 (I£). 
It is not difficult to see that due to the structure of the Einstein equations, Proposition 46 
also holds when all the barred and unbarred quantities are exchanged. The proof is exactly 
analogous to that of Proposition 46 and will be omitted. 

Proposition 47. The following L 2 estimates for the curvature components hold: 

^LlL\S)) 



< 



i<2 



+ 
+ 
+ 
+ 
+ 



LHSo,u)> 

(J]V i (p, ( T))( Yl ^ 1 V 42 (^ + ^)V l3 (p,a))|| i i iiil(5) 

i<2 ii+i 2 +i:j<2 

(£vU*))( E v li rv j ^v i4 (^ + ^))IUiL 4 L 1{ 5) 

i<2 n+42+i3+*4<3 
i<2 ji+i2+i3+*4<2 

(E V ^)( E ^^^(p^MlililHs) 

i<2 i 1 +i 2 +h<2 
i<2 ii+i2+«3<l 

(E V ^ E ^H i2 ^ i3 {trx + ^s)V u (trx + Tp H ))\\ L iL k LHs) 

i<2 ii4-i 2 +i3+H<3 



We now control all the error terms in the energy estimates. Introduce the bootstrap 
assumption: 

K < A 4 . 

First we estimate p and a in L 2 (H U ) and j3 in L 2 (H U ). 



(A5) 
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Proposition 48. There exist €q = eo(Oo, IZo, A 4) such that whenever e < eo, 

J2m i (p^)\\^Lims) + WV'PW WW)) < c(o ,n ) 



i<2 



Proof. We control the six terms in Proposition 46 We first estimate the term (3tpHipH_, i-e., 
the last term in the expression in Proposition 46" As we will see, this is the most difficult 
term because all three factors can only be controlled after taking the L? norm along one of 
the null variables. 

i<2 ii+«2+i3+j4<3 

<e^(J2 I|V^IU-^l»(5))( llV n # 2 V i3 ^V i4 ^|| L|L 2 L2(5) ). 

n+«2+«3+*4<3 



i<2 



Since we have a small constant ea, we only need to bound this by a constant depending on 
A4. The first factor is bounded by A4 by the definition of the norm 1Z and the bootstrap 
assumption (|A5|). We now look at the second factor. 



J2 V^V^V^/dkws) 

)) (Yl 1 1 V * 2 ^ I Il^lIlHs)) I \iPh\ \lIL^L°°(S) 

i 2 <3 

(5) ) I \^H\ |iflL»L<»(S)(5^ 1 1 V l2 ^| | L ^ L 2 L 2 (5) ) 
i 2 <3 

(S ))(I] HV i3 ^IU-L|L4(S))(^ l|V i4 V £ ||L^g=i^ 



<c(£im 


1*1 v 
\L°°LtfL°°(S)< 


ii<3 




+o(£i 




ii<3 




+ 0(1)1 


■0 1 L°° L°° L°° 


ii<3 




+ o(£l 


HM 1 L°° L°° L°° 


ii <1 




+o<Ei 

il<l 





i 2 <2 



i 3 <2 



)) I IV'hI Us°l|l°°(5) I \^h\\lil^l^(s)C^2 1 1 Vi2 ^l |ls°l^l 2 



2 (S), 



i 2 <2 



(J2\\V l2 ij\\ L ~L?LHS)) 
i 2 <l 



+ C| |Vh| U»£|l<»(s)(E 1 1 V ^l kxg> lhs))(Y1 1 1 v * 



|LS°L-L4(S); 



12<1 



By Propositions [36] and [43] we have 

II Yl vH^HhV^hIIlilil^s) < c(o ,n , a 4 ). 

«l+« 2 +«3+«4<3 

Thus 

|| £v*£ £ V^V^V^idk^*^ < C(0 ,^o, A 4 )eJ 

i<2 ii+i 2 +« 3 +j4<3 

We next consider the following four terms from Proposition [46] 

||(£v<(p,a))( £ ^V-(^ + ^)V J3 (p,a))|| LiLiil(s) , 

U+«2+«3<2 



i<2 



i<2 ji+«2+*3+M<3 
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||(2>'£)( E rv i2 VV i3 (p^))iUi ilLl( 5), 

i<2 ii+j 2 +i 3 <2 
i<2 ii+i2+*3<l 

Since ^ satisfies stronger estimates than either or ipjj] and p, a satisfy strong esti- 
mates than either Vipn or Vipn, we can bound these terms exactly the way as above by 
C(O ,n ,A 4 )el 

We are thus left with the term 

||(^V i (p )( r))( Yl ^V^^xV^xIIl^w). 

i<2 ii+i 2 +i3+i4<2 

Since x can only be controlled after taking the L? norm in u, we must bound the curvature 
term V l (p, a) in L 2 (H). Nevertheless, we get a smallness constant in this estimate: 

||(^V i (p ) a))( Yl ^ i2 ^ i3 X^ i4 mhLlLHS) 

i<2 ii+i 2 +i3+i4<2 
i<2 ii+i2+i3+«4<2 

<Ce^(£ £ 1 1 V^| |^ L . i2(5) ) (E 1 1 V i3 %i I W i2( 



2 (S)y 



ii<2i 2 <3 i 3 <2 
1 



<L7(O ,^ )e2A 4 . 
Therefore, 

Ed I V ^> *) I I + W^iWl^Lims)) <K + &C(O , Ko, A 4 ). 

i<2 

Thus the conclusion follows by choosing e to be sufficiently small depending on Oq,TZ and 
A 4 . □ 

We now estimate the remaining components of curvature: 

Proposition 49. There exist €q = eo(Oo, TZq, A 4 ) such that whenever e < eo, 

E(ii v ^H^^ 2 (5) + \\v*(P,*)\\igww) < c(o ,n ) 

i<2 

Proof. In order to prove this estimate, we heavily rely on the bounds that we have already 
derived in Proposition 48 for V l (p, a) and V*/3. In particular, we need to use the fact that 
those estimates are independent of A 4 . In order to effectively distinguish the norms for the 
different components of curvature. We introduce the following notation: 



RuiP] :=Eil V ^lk^, a )> 

i<2 

Ru\p,d-] ■= Ell^^'^H^ 2 ^*)' 



i<2 



Ru[p,o-) := Ell^^'^H^ 2 ^) 



i<2 



NO 
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i<2 



We now proceed to proving the proposition by controlling the six error terms in Proposition 
47 We start with the first, second, fourth and fifth terms: 



J>*(p,a))( rV J2 (^ + ^)V i3 (p,a))|| L i iiLl(5) 

ii+*2+«3<2 



i<2 



and 



(EVfam e v ii ^v^v i4 (^+^))iu i L i L 



i(S) 



i<2 



«l+«2+i3+«4<3 



and 



(S) 



i<2 



and 



;^V il /3)( £ V i2 ^ 3 V i4 ^V i5 trx)|| L i^L 



i(5)- 



ii<2 



«2+*3+«4+«5<3 



In these terms (5 or appears at most once. Therefore, after applying Cauchy-Schwarz in 
u and putting or ipx in L^, there is still an extra smallness constant ez. These terms can 



be estimated in a similar fashion as in Proposition 48 by C(Oq, TZq, A 4 )ea . We then look at 
the last term in Proposition [47} 



E>V?)( E V Jl ^ 2 V 43 (trx + ^)V i4 (tr X + ^))IUiLiL 



i(5)- 



(45) 



i<2 



il+«2+i3+M<3 



We first look at the case where we have ftipHipH- Among these terms, there are two possi- 
bilities: the case where (tr%, ip H ) has at least 2 derivatives and the case where (tr%, ip H ) has 
at most 1 derivative. For the term where (tr%, ipn) has at least 2 derivatives, we have 

U+J2+i3+*4<3,2<i3<3 



< / \\v 2 

'o 



? U ',«)(E 



ii<l 



X 



E ll Vl2 ( tr ^'^)H^ 2 (^,j)(E liV i3 (trx,VH)IU»x4(^, 



2<i 2 <3 



i 3 <l 



< / C(Oo,7^o)(l + 7^ u 4/3] 2 )(^||V ^ (trx,^)||L 5f L 4 (5 ^l ,J)^ , 



by Propositions 36 and 43 For the term where (tr%, ipa) has at most one derivative, notice 
that the estimates for V^(trx, ipn) m L 2 in Proposition 36 depends only on initial data and 
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the bounds for V 3 (tix i[>h) in Proposition 43 depends only on initial data and TZu[(3}. Thus, 

ii(E v ^)( E 



V^V^trx, <MV i4 (trx, ^h))\\ip u ^3) 

i<2 ii+i 2 +j3+i4<3,i3<l 
il<2 i 2 <2i3<3 

x (E l|V l4 (tr X ,^)|| L|L$fL2(s) )(E \\^ i5 (tTX^H)\\ L ^LlLHs)) 

14<2 «5<3 

<C(O ,n )(l + n u [P]){l + Ku[(3}). 



< 



Therefore, ( |45[ ) can be estimated by 

H(E V ^)( E V ll # 2 V J3 (trx + ^)V 44 (tr X + ^))||LiL i L 1( 5) 

C2 ii+i2+i3+«4<3 

C(0 ,Ko)(l + n u ,{{3] 2 )(J2 HV^trx,^)!!^^^, S)du' 
+ C(O 0l K )(l + 7^[/3])(l + TZu[§\). 

We note explicitly that it is important that we do not allow all terms of the type ipH^H^ 
but only allow ipn^H^X since by Proposition 43, 03,2 [trx] can be controlled by a constant 
depending on initial data, but the bound for Os^iVyV] depends on 1Z. 
We are thus only left with the term 

1 1 (E V 41 (p, &))( E V^V^Hh^Hh) \\LlL k L H S), 

h<2 i2+«3+*4+«5+*6<2 



i.e., the third of the six error terms in Proposition |47| Since ipu does not enter with three 
derivatives, it can be estimated using Proposition [36 by 

1 1 (E V 11 (p, a)){ E V* 2 # 3 V^V^V^h) I Ui^(s) 

<(Eiiv ii (p^)iu^ 2 (5))( E 

il<2 j2+«3+*4+i5+«6<2 

<(En vil (p^)iu^(5))(EEii v ^n 

i2<2i 3 <3 



VH^ u ^HhV 1 ^h\\lililhs)) 

xEii v "«WW 



Il-ls°l 2 (s) 



ii<2 



i 4 <2 



Therefore, we have 

ft u [/?] 2 + 7^[p,a] 2 

<C7(0 ,^o)(l + e5C(Oo,^o,A 4 )+ / (7M/?]) 2 (E ^H^L^Jdv!) 



i<l 
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Applying Cauchy-Schwarz on the last two terms and absorbing ^{1Z U [(3] 2 + TZu[p, &} 2 ) to the 
left hand side, we have 

<C(O ,K )(l + e 1 *C(O ,K ,A 4 )+ / (KM) 2 (J2 \\V^ h \\l^ { s u , jW + K U [§} 2 ). 

Jo ^ 



Gronwall's inequality implies 

n^ + n^d-] 2 



<C(O , Ko)(l + e^CiOo, Ko, A 4 ) + TZ^ff ) exp( H V 1 1 V^IUy^, „))<*«') 

Jo i<i ~ 



By Proposition 36 



exp( 



i<l 



By Proposition |48 
Therefore, 



\\V^H\\L~L H s u ,Jdu') < C(G ,Ko). 
llu[(3} 2 <C(O ,7l ). 



Kuffi + KJip,*] 2 < C(O ,K )(l + e5C(0o,fto,A 4 )). 
Taking e sufficiently small depending on Oq, TZq and A 4 , we conclude that 

n u {(3} 2 + 1lu[p,a} 2 < C(O ,K ). 



□ 



Propositions 48 and 49 together imply 
Proposition 50. There exists €q = (Oq,TZq) such that whenever e < eo, 

n<c(o ,n ). 

Proof. Let 

A 4 »C(Oo,^ ), 

where C(Oq,TZq) is taken to be the maximum of the bounds in Propositions 48, and 49 



Hence, the choice of A 4 depends only on Oq and TZq. Thus, by Propositions [48J and|49J the 

on 

□ 



bootstrap assumption (A5) can be improved by choosing e sufficiently small depending on 
Oq and TZq. 



This concludes the proof of Theorem |4} 

7. Nonlinear Interaction of Impulsive Gravitational Waves 

In this section, we return to the special case of the nonlinear interaction of impulsive grav- 
itational waves, thus proving Theorem [T] Recall in that setting we prescribe characteristic 
initial data such that on H (0,y^) (resp. H (0,u*)), x (resp. x) is smooth except on a 
2-sphere S 0>w (resp. S Ust0 ) where it has a jump discontinuity. Thus the curvature in the 
data has delta singularities supported on S 0tVi and S Ust0 . 

Such initial data set can be constructed by solving a system of ODEs in a way similar to 
the construction of the initial data with one gravitational impulsive wave in [25]. Moreover, 
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one can find a sequence of smooth characteristic data that converges to the data for the 
colliding impulsive gravitational waves. We refer the readers to [25] for more details. 

With the given initial data, Theorem [3] implies that a unique spacetime solution (M.,g) 
to the Einstein equations exists in < u < and < u < u^. Moreover, using the a priori 
estimates established in Theorem |4| we can show that the sequence of smooth data described 
above gives rise to a sequence of smooth spacetimes that converge to (Ai,g). 

In this section, we prove that in addition to the a priori estimates proved in Theorem |4j 
the colliding impulsive gravitational spacetime (Ai,g) possesses extra regularity properties 
as described in parts (&),(c) of Theorem [l] We give an outline of the remainder of the 
section: 

Section 7.1 We show the first part of Theorem Tic), i.e., that (3,p,a,(3 can be defined 



in L 2 U L 2 U L? '(S) . This follows directly from the estimates in the proof of Theorem p 



Section |7.2 We prove the second part of Theorem []Jc), showing that the solution is 



smooth away from H_ u U H Ug . 

Section 7.3 We establish Theorem ^b). We define a and a in the colliding impulsive 



gravitational spacetime and show that they are measures with singular atoms supported on 
H u and H Us respectively. This shows that the singularities indeed propagate along the null 
hypersurfaces H Us and II u . 

7.1. Control of the Regular Curvature Components. 

Proposition 51. All the curvature components except a and a are in L 2 U L 2 U I? [S) . 

Proof. It follows directly from the proof of Theorem[3 that (3, p,a, (3 £ L\L\L?[S\ It remains 
to show that p, a are in L 2 Ju 2 u L 2 {S). Recalling the definition of p and a, it suffices to show 
that XX is i n L 2 a L 2 l L 2 (S). This follows from 

WxxWlilil^s) < \ \x\\l~lilhs)\\x\\lil^lhs)- 

□ 

7.2. Smoothness of Spacetime away from the Two Singular Hypersurfaces. In 

this subsection, we prove that in the case of two colliding impulsive gravitational waves, 
the spacetime is smooth away from the null hypersurfaces H Us and B_ u . For u < u s or 
u < u s , this follows from the result of We will therefore only prove the statement for 
{u > u s } H {u > u s }. 

Proposition 52. The unique solution to the vacuum Einstein equations for the initial data 
of colliding impulsive gravitational is smooth in {u > u s } C\ {u > u s } 

Proof. We establish estimates for all derivatives of all the Ricci coefficients. We prove by 

induction on j, k that V'V^V*^, ^h, iPh) is in L^L^L 2 (S) for all i, j,k. 

1. Base case: j = k = 

1(a). Estimates for ip, p and a 

For i < 2, V l, is in L^L^L 2 {S) by Theorem |4 Using exactly the same arguments 
but allowing more angular derivatives in the initial data, it is easy to show that V l, is in 
L™L™L 2 (S) for all i. 

Similarly, an adaptation of the arguments in Theorem [4] imply that V 4 /5 and V*cx are in 
L™L™L 2 (S) for all i. 
1(b). Estimates for tpn and tpn_ 
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The a priori estimates given by Theorem [4] only imply that for % < 2, 

Vtyir G LlL™L 2 (S) and Vtyg G L 2 u L™L 2 {S). (46) 

Applying a simple modification to the proof of Theorem [4] with more angular derivatives in 
the initial data, it is easy to show that (46) holds for all % > 0. In order to improve this to 
a bound in L 2 (S) , we need to use the fact that we are away from the hypersurfaces 

H Us and H_u s - 

We first prove estimates for \. Consider the equation 

V 3 X + ^trxx = V®t? + 2ux - ^tr^X + ?7®?7- (47) 

Since we know that the initial data on H n {u > u s } are smooth, V l x is in L^L 2 (S 0t u). 
Using the control that has already been obtained, we integrate (47) to show that V*x is in 
L™L™L 2 (S) for u > u s for all i/ 
Similarly, using 

V 4 X + ^XX = V% + 2ux - ^trxx + 77%, 

we show that V*x is in L^°L2°L 2 (5') for u> u s for all i. 

The estimates for p and <r together with the bounds for x an d X imply that V l p and W 
are in L^L^L 2 (S). Now integrate 



and 



3 1 11 

V 3 w = 2uu+ -\v-v\ 2 + -(v-v) -(V + V)- g\v + v\ 2 + 



3 1 11 

V 4 w = 2ww + -\r] - r/| 2 - -(77 - 77) • (77 + 77) - -\rj + rf + -p, 



to show that V % u and VV are in L^L^L 2 {S) for all i. 
2. Induction Step 

We now proceed to the induction step. Assume V^VgV^V*, i>H, ^h) is in L 2 (S) for 

all i, for all j <J and for all k < K. We show that V'V^+M^, Vte, ^jr) is in L™L™L 2 (S) 
for all % and for all k < K. A similar argument then shows that +1 (ip, ipHi iI>h) is 

in L^°L^ C 'L 2 (S') for all z and for all j < J. This completes the induction step and prove the 
proposition. 

For ipu and ip that satisfies a equation, we commute the equations J times with V3, 
k < K times with V4 and differentiate i times with V to get 

V i V 3 J+1 V 4 fc (V,^) = - 
where every term on the right hand side can be controlled in L 2 L 2 L 2 (5') by the induction 
hypothesis. Thus we can estimate directly to show that V l V^ +1 V^(tp, ipa) is in L^L^L 2 (S) 
for all i and all k < K. 

For ipH_ and i[) that satisfies a equation, we commute the equations i times with V, 
J + 1 times with V3 and k < K times with V4. Then we have 

V 4 (V i V 3 J+1 V^^,^)) = (^,^,^)v l v 3 J+1 v^(^,^) + ... 

where ... are terms that can be bounded in L 2 L 2 L 2 (S) using the induction hypothesis. Thus, 
by GronwalPs inequality, V i V 3 7+1 Vt(V' ) V'g) is in L u L u L2 (S) for all i and all k < K. □ 
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7.3. Propagation of Singularities. We first show that a and a can be defined as measures. 
We take the null structure equations 

a = -V 4 X - tixx ~ 2^X, (48) 

QL = - V 3 v - tixx ~ 2wx 

as the definitions of a and a. In view of the fact x an d X are n °t differentiable, a and 
a cannot be defined as functions. Nevertheless, we will show that they can be defined as 
measures. By (48), if a is smooth, we have 

f- d 

a(u,u,'d)du = J (ti^ 1 — x + tixx + 2ux)(u,u,'d)du. 
Integrating by parts and using 

d 

—Q- 1 = 2u, 
ou 

we derive 

a(u,yf,-d)dv! = (n _:L x)(«, u,$) - (fi -1 x)(it,u = 0,0) + / (tr xx)(u,u',$)dyf. 



For every u ^ u s , the right hand side is well-defined. For each u, i? G § 2 , we define a as a 
measure such that 

«([0,w)) = (tt~~ 1 X)(u,%'$) - (ft -1 x)(ii,u = 0,t?) + / (trxx) (u,u,$)du for u ^ u s . 

Jo 

By continuity, we have 

a([0,u s )) = lim (n^xXu,^) - (fi,~ 1 X )(u, u = 0,$) + (tr X x)(u, yf^)dvf. 



a([0,u)) = (r^KM,!?) - (Q-'x^u = 0,u,$) + / (fi-V Jjx + trxx)^',"^)^', 



This defines a as a measure. Similarly, for each u, $ G S 2 , we define a to be a measure by 

dd A - 

for u ^ u s . By continuity 

«([0,iO) 

= lim (fi- 1 *)^!?) - (n^xXu = 0,«,i?) + r\n- 1 b A J^x + trxxKu',u,$)du'. 

u^u 3 Jo OU 

Remark 1. //we tofce a sequence of smooth initial data converging to the data for nonlinearly 
interacting impulsive gravitational waves, it can be shown that in the spacetimes (Ai n ,g n ) 
arising from these data are smooth and oc n —> a, a n — >■ a weakly, where a and a are as defined 
above. We refer the readers to [25J for details in the case of one impulsive gravitational wave. 

Proposition 53. x ^ s discontinuous across u = u s . Similarly, x is discontinuous across 
u = u s . 

Proof. We focus on the proof for x- The proof for x is similar. Consider the equation. 

V 3 V + ^trxx = V®n + 2u_x - ^trxx + 
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Rewriting in coordinates, we have 

d d 1 ~ 1 ~ 

Consider the coordinate system (u, u, O 1 ^ 2 ) such that 

u = u- b A 6 A . 

Since < Ce, this is well-defined in a local coordinate patch. For fixed u, 8, define Uq to 
be the w value such that (w , u, 0) e H . Thus 

d 1^1^ 
fil (^ + 7 _1 X)X + otrxx = V<g>?7 + 2ur£ - -tr%% + (49) 
cm — 2 — 2 — 

For the initial data, x(«0) is smooth for u ^ u s and has a jump discontinuity for 



u = u s . The conclusion thus follows from integrating (49) and using already established 



estimates. □ 

Finally, we show that a (resp. a) has a delta singularity on the incoming null hypersurface 
H u (resp. outgoing hypersurface H Us ). 

Proposition 54. a can be decomposed as 

a = d~(u s )a s + «r, 

where S(u s ) is the scalar delta function supported on the null hypersurface H_ u , a s = 
a s (ii,$) 7^ belongs to L 2 L 2 (S) and a r belongs to L 2 L 2 L 2 (S). 
Similarly, a can be decomposed as 

a = d~(u s )a s + a r , 

where S(u s ) is the scalar delta function supported on the null hypersurface H Ug , a s = a s {u, ^ 
belongs to L^L 2 (S') and a r belongs to L 2 L 2 L 2 (S). 

Proof. We prove the proposition for a. The statement for a can be proved in a similar 
fashion. Define 



a s (u, $) := lim x(w,w, — lim x(w, w, 

and 

a r := a — S(u s )a s . 

We now show that a s and a r have the desired property. By Theorem |4| a s belongs to 
L 2 1 L 2 (S'). That a s ^ follows from the fact that x has a jump discontinuity across u = u s , 



which is proved in Proposition 53 
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It remains to show that a r belongs to L 2 L 2 L 2 '(S) . To show this, we consider the measure 
of the half open interval [0,w) using the measure a r (u, i?): 

(a r (u,tf))([0,M)) 

= (fi- 1 x)(M,M,^) - lim (ft- 1 ^,^) + lim (fi -1 x)(u, u, 0) - (JT^Xu, u = 0, #) 







<9w 



lim 



(r 1 x)(u,u' ) )?)«'+ lim 



By Proposition 

Thus (a r (u, $)) 
as desired. 



52 



3 



^(m, m, is hi L 2 L 2 L 2 (S) away from the the hypersurface H u . 
~[0,m)) can be expressed as an integral over [0, u) of an L^L^L 2 (S) function, 

□ 



8. Formation of Trapped Surfaces 

We also apply the existence and uniqueness result in Theorem [3] to the problem the 
formation of trapped surfaces. In [8], Christodoulou proved that trapped surfaces can form 
in evolution. This was later simplified and generalized by Klainerman and Rodnianski [19J, 
[20]. These are also the first large data results for the long time dynamics of the Einstein 
equations without symmetry assumptions. 

In all the previous works, the setting is a characteristic initial value problem such that 
the data on the incoming null hypersurface are that of Minkowski spacetime. The data 
on the outgoing null hypersurface, termed a "short pulse" by Christodoulou, are large, but 
are only prescribed on a region with a short characteristic length. The large data on the 
outgoing hypersurface and the small data on the incoming hypersurface together give rise to 
a hierarchy of large and small quantities, which was shown to be propagated by the evolution 
equations. 

In particular, in order to guarantee the formation of a trapped surface, the initial norm of 
X is large on H , and is of size 

_ i 

where e is the short characteristic length in the u direction. Moreover, a has initial norm of 
size 

1 1 1 1 -- 
IMIl2°l°°(S) ~ e 2 - 

It was precisely to offset the largeness of x an d a (and their derivatives) that the data on 
if o were required to be small. 

However, when viewed in the weaker topology L 2 L°°(S'), the initial size for x in [8] is 
bounded by a constant independent of e: 

HxlU|i°°(s) ~ i- 

Our main existence result applies for initial data such that x an d its angular derivatives are 
only in L 2 U L°°(S) without any requiring any smallness for the data on H^. In particular, 
no assumptions on a and its derivatives are imposed. Using this theorem, we obtain the 
following extension to the theorem in [8] , [IE] : 
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Theorem 5. Suppose the characteristic initial data are smooth on H_ Q for < u < u* and 
satisfy the following inequality: 




J exp(-y trx(u'\$)du")(-2K + 2div ( + 2\C\ 2 )(u' ,$)du' 



(50) 



trx{u', d)du')trx{u = 0, $) 



for every $ G S> 2 . Then there exists an open set of smooth initial data on H such that the 
initial data do not contain a trapped surface while a trapped surface is formed in evolution. 

More precisely, there exists e > sufficiently small such that if the characteristic initial 
data on H D {0 < u < e} are smooth and satisfy 

52\\V*x\\i&>(S) < C (51) 
and the following two inequalities are verified for every d G § 2 , 

\x\ 2 Wu 



1 f u 

>exp(- / trx(u,u = 0,$)du) ( 52 ) 



trx(u 



0,m = 0,?9) + ^ exp (\f Q tr xdu")(-2K + 2div C, + 2\C,\ 2 )du\ , 



and 

\x\ 2 ($)du < trx(u = 0,u = 0,??), (53) 



then there exists a unique spacetime (A4,g) endowed with a double null foliation u, u and 
solves the characteristic initial value problem for the vacuum Einstein equations in the region 
< u < u*, 0<w<e. Moreover, H fl {0 < u < e} does not contain a trapped surface and 
S Ut ,t is a trapped surface. 



Remark 2. (50) holds in particular on a regular null cone, such that 

2 2 

\\trx- -, tr X + -,(, VC,pIU~l°°(s„,o) < C, 

for u* sufficiently close to the vertex. 

In particular, this implies the celebrated theorem of Christodouloi 

Corollary 55 (Christodoulou). // the characteristic initial data on if is that of the trun- 
cated backward light cone 

{u = t + r = 0, 0<u<l 



in Minkowski space, then for e sufficiently small, if the data on Hq satisfy (51), (52) and 
(53), then there exists a unique spacetime (A4,g) endowed with a double null foliation u, u 



8 The original theorem of Christodoulou in [8] constructs a spacetime from past null infinity. Here, we 
retrieve only the theorem in a finite region. Nevertheless, the infinite problem can be treated as in [5] once 
the finite problem is understood (see also a simplified treatment in [3]). 

9 Here, we adapt the notation that u = t — r + 2, u = t + r. Therefore, < u < 1 corresponds to the 
i-range —2 < t < —1. 
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and solves the characteristic initial value problem for the vacuum Einstein equations in the 
region 0<u<l,0<u<e. Moreover, H R {0 < u < e} does not contain a trapped surface 
and Si >e is a trapped surface. 

We now begin the proof of Theorem [5] We need the following series of Propositions. First, 
it is easy to see using the null structure equations and Bianchi equations on H that the 
assumptions for Theorem [3] are satisfied. 

Proposition 56. Given the assumptions for Theorem^ the initial data satisfy the assump- 
tions of Theorem^ Therefore, using the conclusion of Theorem^ a unique spacetime (M.,g) 
endowed with a double null foliation u, u and solves the characteristic initial value problem 
for the vacuum Einstein equations in the region 0<u<u*,0<u<e. Moreover, all the 
estimates in Theorem hold. 



Proof. Since the initial data on H_ is smooth, there exists c and C such that 



j<3 



E (||VV||z^(s„, o) + ||VtyldlL»CSo)) < C ' 



i<3 



E ll^ll^ca) + E 11^*11^(^0) < a 

By (51), x satisfies the bounds in the assumptions of Theorem [3] By the null structure 
equations and the Bianchi equations, for e sufficiently, all the norms for the initial data on 
H Q in the assumptions of Theorem [3] are controlled by a constant independent of e. □ 



We now use the a priori estimates derived in Theorem |4j together with (52) and (53) to 
show that the initial data do not contain a trapped surface and that a trapped surface is 
formed dynamically. We first show that there are no trapped surfaces on H : 

Proposition 57. There exists e sufficiently small such that for all d, 

trxiu = 0, u, i?) > for all u G [0, e] . 

Proof. On Hq, since Q = 1, tr% satisfies the equation 

^tr X = -^(tr X ) 2 -|x| 2 . 



u i 

trx(w = 0,u,&) = trx{u = 0,u = 0,i?) - / (-(trx) 2 + |x| 2 )(«', fydy!. 



Integrating the equation for tr%, we have 

= n. - I ( 

'2 
Hence 

trx(u = 0, u, •&) > trx(u = 0, u = 0, i9) — / |x| 2 (w', ~d)du — Ce. 

Jo 



(53) implies that for every i?, 

tr X (u = 0,u = 0,t?) > / \x?{y!^)dv!. 



o 
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Therefore, for e sufficiently small, 

tvx(u = 0,«,i?) > 0, 

for all u G [0, e] □ 

We now prove that S u , >e is a trapped surface. As a first step, we solve for trx on S Uj o on 
the initial hypersurface H^. 

Proposition 58. On the initial hypersurface H_ Q , trx{u,u = 0, $) is given by 
trx(u, u = 0, #) 

= exp(— - / trxdu')(trx(u = 0, u = 0, $) 
2 Jo ~~ 

exp(- / ^?/)(-2ir + 2dwC + 2|C| 2 )c/n / ). 



+ 



Proof. On if , since f2 = 1, we have 



^-trx + ^tr X tr X = 2p + 2div C + 2|C| 2 . 
aw 2 — 



Substituting the Gauss equation 



K = -p- ^tr^trx, 



we have 

|-tr X + ^trxtrx = -2K + 2div C + 2|C| 2 . 
The conclusion follows easily. □ 
We compare J™* trx(u', u, $)du' and J Q U * trx(u', u — 0, $)du' in the following Proposition: 
Proposition 59. For every u G [0, e], we have 

1/ trx(u',u,$)du - trx(u,u = 0,$)du\<Cek 
Jo Jo 

Proof. Consider the equation 

V^rx = — -tr%tr% + 2utix + 2p + 2div rj + 2\r]\ 2 . 

Writing V4 = and integrating, it is easy to see that by the estimates in Theorem 

we have 

|tr%(w, u, -d)dv! — trx(u,u = 0, $)| < Ce^ for all u for all $ G § 2 . 
This implies the Proposition after integration in u. □ 

Using Proposition 59, we compute f Q Ixl^du for every u and every ■& G § 2 : 

Proposition 60. For every -d in S Uje , the integral of\x\~ along the integral curve of L through 
(u, 1?) satisfies 

f u f € 1 

|x| 2 (w, u, d)du > exp(— / tr~x(u' ,u = 0,$)du') / |x| 2 ( M = 0, u, d)du — Ce? . 
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Proof. Fix '&. Consider the null structure equation 



V 3 X + ^trxx = V<g>r? + 2ux - ^trxx + V®V- 



Contracting this two tensor with \ using the metric, we have 
2 



19 9 1 1 

+ &A ^7T4 + o tr xlA^ - 2 ^lxl7 = X(^®V ~ ^XX + V®v)- 



In coordinates, we have 

-(- + b A - 
Using 

u = --V 3 (logfl), 

we get 



^0-21-12 

A I <y 



f2 2 exp(— / Vttixdu)— ( exp( / Qtrxdu')£l 
d dVt 1 



(54) 



Let 



/« 9 <9fi 1 

mrxdttOC-^^jIxl? - + 2fix(V®r/ - 2fi-tr X X + 2^77)). 

Using the equation 

cm 

the estimates for Q and £ and the fact that b A = on iJ , we have a uniform upper bound 
for b A : 

\b A \ < Ce. 

Thus, 

\\F\\lilIl°°(S) < Ce*. 



Therefore, by (|54|), we have 

exp( / Q(u',u)trx(u ', u)du'){l~ 2 (u,u)\x\^,(u, u) 



>lxl> = 0,u)-C||^(w)lk£o 



The proof of Proposition [T] implies that 

||^ _ 1\\l°°l°°l°°(S) < Ce 2 . 

This, together with Proposition |59j gives 

Q~ 2 (u, u) exp( J Q " Q(u', u)trx(u', u)du') 
exp( J" m = 0)trx(w', w = 0)ofat') 



II < Ce 2 . 



92 



JONATHAN LUK AND IGOR RODNIANSKI 



Therefore, 

exp( / trx(u',u = 0)du')\x\ 2 (u, u) 
Jo 

>lxl?(« = 0,w) - Cel\x\*(u,u) - C\\F(u)\\ LlL oo {s) . 
Taking the L\ norm, we get 

exp( / trx(V, u — 0)du') / \x\j(u,u)du 
Jo ~ Jo 

> \x\ 2 1 (u = 0,u)du-Ce 1 i \x\y(u,u)du- C\\F(u)\\ L ^ L i L oo {s) 
Jo Jo 

> [ \x\ 2 1 {u = 0,u)du-Ce*. 
Jo 

This allows us to conclude the formation of trapped surface: 



□ 



Proposition 61. Given the assumptions of Theorem^ for e sufficiently small, trx < 
pointwise on S Utte . In other words, S Uit , e is a trapped surface. 



Proof. By Proposition 58, we have 



tix(u„,u = 0,t?) 



ex P(-^ / tTxdu')(trx(u = 0,u = 0,1?) 

2 Jo _ (55) 



+ 



pu* i pu' 

/ exp(- / trx^")(-2K + 2divC + 2|C| 2 )). 
Jo 1 Jo 



By Proposition 60 



\x\y(u*,u,$)du> exjp(— / trx(u , u = 0, d)du') / |x| 2 (w = 0, u, fi)du — Ce 2 . (56) 
o Jo ~ Jo 

Using the equation 

V 4 trx = -^(trx) 2 - \x? - 2wtrx, 

which can be written in coordinates as 

_ d 1 

we have 

trx(«*,w = e, i?) < trx(u*,U = 0,1?) - / \x\ 2 (u if ,u,^)du + Ce^ . 

Jo 
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Therefore, using (55) and (56), we have 



< 



trx(«*,M = e,i?) 
1 m, 



tTx(u',u = 0,$)du') 



i p 



x I trx(u = 0,u = 0,1?) + / expQ I trxdu")(-2K + 2div ( + 2\(\ 2 )du' 



exp(- / trx(w',i| = 0,i?)di/) / \x\l(u = 0,u,'&)du + Ce* . 



Since by (52), for all i?, 



trx(w = 0,u = 0,i?) + J exp(-y trxdw")(-2^ + 2div ( + 2 [CI 

<exp(— - / trx(u, u = 0, $)du) f \x\*(u = 0, u, $)du, 
2 Jo ~ Jo 

e can be chosen sufficiently small so that 

trx(u*,u = e, i?) < for every i?. 



2 )dv! 



□ 
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